MZ-TH/00-52 
DESY-/01-065 
hep-th/0012189 
December 2000 



New results for the e-expansion of certain 
one-, two- and three- loop Feynman diagrams 



A. I. Davydychev^'l] and M. Yu. Kalmykov^-Q 



"'Department of Physics, University of Mainz, 
Standing erw eg 7, D-55099 Mainz, Germany 

^ DESY-Zeuthen, Theory Group, Platanenallee 6, 
D- 15738 Zeuthen, Germany 



Abstract 

For certain dimensionally-regulated one-, two- and three-loop diagrams, problems 
of constructing the e-expansion and the analytic continuation of the results are studied. 
In some examples, an arbitrary term of the e-expansion can be calculated. For more 
complicated cases, only a few higher terms in e are obtained. Apart from the one- 
loop two- and three-point diagrams, the examples include two-loop (mainly on-shell) 
propagator- type diagrams and three- loop vacuum diagrams. As a by-product, some 
new relations involving Clausen function, generalized log-sine integrals and certain 
Euler-Zagier sums are established, and some useful results for the hypergeometric 
functions of argument 4 are presented. 



*0n leave from Institute for Nuclear Physics, Moscow State University, 119899, Moscow, Russia. Email 
address: davyd(Q)thep. physik.uni-mainz.de 

tOn leave from BLTP, JINR, 141980 Dubna, Russia. Email address: kalmykov@ifh.de 



1 Introduction 



One of the most powerful tools used in loop calculations is dimensional regularization [|l|. In 
some cases, one can derive results valid for an arbitrary space-time dimension n = A — 2e, 
usually in terms of various hypergeometric functions. For practical purposes the coefficients 
of the expansion in e are important. In particular, in multiloop calculations higher terms of 
the ^-expansion of one- and two-loop functions are needed, since one can get contributions 
where these functions are multiplied by poles in e. Such poles may appear not only due to 
factorizable loops, but also as a result of application of the integration by parts or other 
techniques [^. 

One of the important examples are massless propagator-type diagrams occurring in the 
renormalization group calculations. By now, the structure of the terms of the e-expansion 
of such diagrams is well understood. As a rule, the occurring transcendental numbers can 
be expressed in terms of multiple Euler-Zagier sums 



k 



•^7 

ni>n2>...>rafc>0 j=l "'J 



where Oj = ±1 and Sj > 0. For lower cases, these sums correspond to ordinary (^-functions. 
For studying the higher-order sums analytical and numerical methods have been recently 
developed It was demonstrated that there is a "link" between the quantum field 
theory and the knot theory [|: some Feynman diagrams can be connected with knots, so that 
the values (Euler-Zagier sums) of Feynman diagrams are also associated with knots. 

In multiloop massive calculations in QED and QCD [0 within the on-shell scheme, new 
constants appear, which are related to polylogarithms. A detailed description of the basis 
of this type is presented in Ref. Two- loop vacuum diagrams with equal masses f^, |10 



yield the transcendental number CI2 (^|), where Clj (9) is the Clausen function ( |A.4|) . Some 



useful properties of this function are collected in Appendix A (see also in [TT[). The same 
constant CI2 j appears in the one-loop off-shell three-point diagram with massless internal 
lines, in the symmetric case when all external momenta squared are equal, see in |12 



To classify new constants appearing in single-scale massive diagrams, Broadhurst has 



introduced in Ref. 13| the "sixth root of unity" basis connected with 





E (1-2) 



^7 

ni>n2>...>nj,>0 j=l "'J 



where A = exp (^i|^ and pj G {0,1,2,3,4,5}. For particular cases pj G {0,3} it coincides 
with the Euler-Zagier sums ( |1.1|) . For both definitions ( |1.1|) and (|1.2|) , the weight can be 
defined as J2i=i Si, whereas the value of k can be associated with the depth (see in ||F 



'^An on-line calculator for the Euler sums, with an accuracy of 100 decimals, is available at 



http : //www. cecm. sf u. ca/projects/EZFace/ . 

•^A collection of results related to this connection can be found on Dirk Kreimer's home page, 



http : //yan. open. ac . uk/ ~dbroadhu/knf t . html 



http : / / dipmza . physik . uni-mainz . de/ kreimer/ps/knf t .html or on David Broadhurst's home page 
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One of the remarkable results of Ref. |T^ is that all finite parts of three-loop vacuum 
integrals without sub divergences, with an arbitrary distribution of massive and massless 



lines, can be expressed in terms of four weight-4 constants: ^4, CI2 ) 



f/3,1 = -2Li4 
(see Eq. (55) of |jl3|) and, finally. 



+ iC4-3^1n'2 + iC2ln'2 



p>fc>0 



cos 



(1.3) 



:i.4) 



which is an essentially new constant. The same constants appear in the (2 — 2£:)-dimensional 
Feynman integrals [Q, as one should expect due to algebraic relations between diagrams 

with shifted dimension 0. We also note that CI2 (|) appears in the two-loop non-planar 
three-point diagram [|l^], when internal lines are massless, whereas all external momenta 
squared are off shell and equal. Unfortunately, a large number of elements (more than 4000) 
makes it difficult to define the complete basis of "sixth root of unity" at the weight 4. So 
far, only the cases with the depth k <2 have been examined [113 . 



In Refs. [^, , the central binomial sums were considered, 

V 1 



S{a) ^ E 



:i.5) 



It was shown that these sums are connected with the multiple Clausen values, which are 



defined as the real or imaginary part of the mult i- dimensional polylogarithm |18 



Li 



11, 



.,a,(^)= E 

"l>"2>--->"fe>0 



k 

n 



ai 
3=1 



n 



:i.6) 



taken aX z = exp fi|j (i.e., the "sixth root of unity" )0. In particular (see Theorem 1 in |]T 



Cs - t|C5,3 - if CsCs + iC2C; 



3' 



:i.7) 



5(8) = -47rIm[Li6,i(e-/=^)j , ,33,3, 

where we use the following short-hand notation (see Eq. ( |1 . 1| ) ) 

=C(si,...,Sfc;l,...,l) . (1.8) 

The constant (^5^3 has occurred in the 6-loop calculation of the /3-function in 0^-theory [jl9| 
(this was recently confirmed in 0]). At the 7- loop order, a new transcendental number 
arises, ^3,5, 3 UTQ). The ^5^3 and another constant, ^7,3, appear in the calculation of anomalous 
dimensions at 0{1/N^) in the large- iV limit A remarkable property of these constants is 

•^An on-line calculator for the multiple Clausen values, with an accuracy of 100 decimals, can be found 



at tittp : //www. cecm. sf u. ca/projects/ezf ace+/ 



3 



their connection with knotsQ |1^, namely, the torus knots 819 and IO124 are associated 
with [29^8 ~ 12C5,3] and [94^10 — 793^7,3], respectively, whereas ^3,5,3 is associated with a 
certain hyperbolic knot [^, |2l[| . 

To predict types of functions (and the values of their arguments) which may appear in 
higher orders of the e-expansion, a geometrical approach happens to be very useful. 
Using this approach, the results for all terms of the e-expansion have been obtained for the 
one-loop two-point function with arbitrary masses p3| , . Moreover, all terms have been 
also obtained for the e-expansion of one-loop three-point integrals with massless internal 
lines and arbitary (off-shell) external momenta and two-loop vacuum diagrams with arbitrary 
masses [0, |2^, which are related to each other, due to the magic connection ||2^. All these 



results have been represented in terms of the log-sine integrals (see in |[TT| and Appendix A.l 
of this paper), whose angular arguments have a rather transparent geometrical interpretation 
(angles of certain triangles). In more complicated cases, like, e.g., the three-point function 
with general values of the momenta and masses, an arbitary term of the e-expansion can be 
represented in terms of one-fold angular integrals whose parameters can be related to the 
angles accociated with a four-dimensional simplex. 

In Ref. [^, the on-shell values of two- loop massive propagator-type integrals have been 
studied, and it was observed that the finite (as e — > 0) parts of all such integrals without 
sub divergences can be expressed in terms of three weight-3 constants, two for the real part, 
(5 and 7rCl2 (f), and one for the imaginary part, 7r(2- 

Furthermore, in Ref. an ansatz was elaborated for constructing the "irrationalities" 
occurring in the ^-expansion of single-scale diagrams involving cut(s) with two massive par- 
ticles. This construction is closely related to the geometrically-inspired all-order e-expansion 
of the one- loop propagator-type diagrams [^, Q , which was also useful when fixing the nor- 
malization factor The procedure of constructing the ansatz is as follows: for each given 
weight j the set {bj} of the basic transcendental numbers contains (i) all products of the 
lower-weight elements {bj^kbk}, k = 1,2, ■ ■ ■ , j — 1 and (ii) a set of new (non-factorizable) 
elements {bj}, which are associated with the quantities arising in the real and imaginary 



parts of the polylogarithms Lij ^e'^^ and Lij (l — e'^^ , with 6* = | or 6* 



2n 
3 



The real and imaginary parts of such polylogarithms can be expressed in terms of the 
Clausen function Clj (6*) (|A.4|) , log-sine integrals Lsj(6') (|A.5|) and generalized log-sine inte- 
grals Lsf\e) Q (see also in ||ll|). Note that Ls2(6') = CI2 (0). The relevant relations are 
collected in Appendix A.l, Eqs. ( [A.2|) and ( |A.7|) . Therefore, in our case the non-factorizable 
part of the basis can be expressed in terms of the functions (|A.4|) , ( |A.5|) and ( |A.6D of two 
possible angles, 6^ = 1 and 9 = It should be noted that this basis is not uniquely 
defined, since there are several relations between polylogarithmic functions Clj (6), Ls j{6) 
and Ls j'\6) of these arguments (see, e.g., Eqs. ( |A.9| )-( [A.11| ) and ( |A.14D of this paper). Af- 
ter excluding all linearly-dependent terms, the basis contains the following non-factorizable 
constants: Lsj(^^^, where j = 3,4,5; Lsj(^|) for j = 2,4,5; Lsf^(^^^ for j = 4,5; and 

Lsf^(^^^. This set of elements will be called the odd basis The numerical values of 

Below we shall discuss connection of these constants with the odd basis. 
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these constants are given in Appendix A of . 

One can use the PSLQ algorithm to search for a hnear relation between the given 
term of the e-expansion of the diagram of interest and the set of numbers {bj}. Using this 
procedure, it was possible to find results for several two- and three-loop single-scale diagrams 
p8[] , some of them have been also calculated analytically in |2^, The constructed odd 
basis has an interesting property: the number Nj of the basic irrational constants of a weight 
j satisfies a simple "empirical" relation Nj = 2\ which has been checked up to weight 4. We 
note that the constant Vs^i (given in Eq. (|1.4|) ) can be expressed in terms of the weight-4 
elements of the odd basis, 

Vs, = i [CI. (f)]^ - i.Ls3(f ) + iCslnS - fiC4 + iLs«(f ) . 

(see Eq. (12) in Ref. p8|]). For the weight 5, 32 linearly-independent elements were elabo- 
rated in Ref. pHf. It should be noted, however, that some of these irrational numbers, like 



Lss (^j; SO far have not appeared in the results for Feynman integrals^. So, the question 
about the completeness of the set of weight-5 basis elements is still open. We are going to 
re-analyze this basis below. 

By analogy with the odd basis introduced in p8|, it is possible to consider the even basis 
connected with the angles | and vr. Apart from the well-known elements tt. In 2, (j and the 

Catalan's constantQG, this basis also contains (up to the weight 5) Lij (^^^ {j = 4, 5, see also 

in §), Lsj(f) (j = 3,4,5)1], CI4 (f) and Lsf^(f). Instead of Li4(|) and Li5(|), one could 

take, e.g., Ls^l^ (^^^ and Ls^^^ (^^^ (or Ls4^''(7r) and Ls5"'^^(7r)), using the relations between these 
elements presented in Appendix A, Eq. ( |A.14 ). For the even basis, the same "empirical" 
relation Nj = 2^ is valid up to weight 4. However, at the weight 5 only 30 independent 
elements have been found so far. This issue will be also discussed below. 

Later, in Ref. [Q, an interesting connection between the function S{a; z) associated with 
the central binomial sums ( |1.5| ) and the generalized log-sine integrals was established, 

^(«; - E = - E , (In zT-'-^ Ls«, (2 arcsin ^] , (1.9) 

where a > 2. In particular, the sums ( |1.5| ) can be represented as 

S(a).S(<.;l)^||gi^-t!lr!Ls...(|). (1.10) 

whereas S{a] 3) is related to Ls^^-* (plus a combination of lower terms with In 3). There- 
fore, S{a\ 1) and S{a] 3) are connected with the terms of the odd basis. Furthermore, consid- 
ering S{a] 2) and S{a] 4) we obtain Ls^^^(^|^ and Ls^^^(7r) (plus a combination of lower terms 
with In 2), respectively, which are related to the even basis. 

^The only known integral containing Ls^^^ {'^) is D5(l, 1, 1, 1, 1,0) (2|. The integral D5(l, 1, 1, 1, 1, 1) is 

reduced to it by using recurrence relations I t is natural to expect that its e-term may contain Lsg^"* (^) • 
•^Note that G = CI2 (f ) = Ls2(§), see in flj. 

^An example of a physical calculation where the constant Ls3(-|) arises is given in [^ . 
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Moreover, in Ref. a more general case of multiple binomial sums has been examined, 

n=l " 

(1.11) 

where Sa{n) = J2]=ij is the harmonic sumQ. In [^, mainly the case k = 1 (related to the 
odd basis) has been treated. A PSLQ-based analysis has shown that not all of sums (|1 . 1 1| ) 
are separately expressible in terms of the basis elements. 

This paper is organized as follows. In Section 2 we examine in detail the ^-expansion of 
the two-point function with arbitrary masses, including the problem of analytic continuation 
of the corresponding functions. In Section 3 we consider some examples of the e-expansion 
of one-loop three-point functions, and study which functions and transcendental constants 
may occur in the cases considered. In Section 4 we consider some physically relevant two- 
and three-loop diagrams. In Section 5 we discuss the results obtained in this paper. There 
are also two appendices containing some further technical details and useful formulae. In 
Appendix A we collect relevant results for the polylogarithms and associated functions. 
In particular, we discuss the relation of the Euler-Zagier sums and the generalized log-sine 
functions. In Appendix B we discuss the expansion of the occurring hypergeometric functions 
with respect to their parameters. 



2 One-loop two-point function 

2.1 General case 

Consider the one-loop two-point function with the external momentum k and masses mi 
and 7712, 



n;i/i,z/2j 



2.1 



[g2 — mfl"^ [{k — qY — mlY'^ 
Using, for instance, a geometrical approach |^2|, one can obtain the following result (for 
unit powers of the propagators, z/i = z/2 = 1): 



J(2)(4-2£;1,1) = in'-'Tie) — \{e + ml-ml)m^'',F, 



1, e 

3 
2 



2 



COS r, 



01 



+ {k^ -ml + 7712)7712 ^ 2F1 



1, E 
3 
2 



2 



cos T, 



02 



(2.2) 



where 2-^1 is the Gauss hypergeometric function. The angles Tq- are defined (see in [^) via^ 



COST, 



k"^ + m{ — 7nl 



01 



2mivk^ 



COSTf 



k'^ — 7nl + 77il 



02 



2m')\/k'^ 



(2.3) 



^Another notation for Sa(n), used by mathematicians, is In particular, for a = 1, Si{n) — Hn- 

When there are no sums of the type Sa{n — 1) or Sb{2n — 1) on the r.h.s. of Eq. (1.11), we shall put a "— " 
sign instead of the indices (a,«) or {b,j) of S, respectively. If the argument (k) is omitted in (1.11), this 
means that the case fc = 1 is understood. 

^ They are related to the angles roi used in Ref. as T-Qi = f — tqi . 



6 



In particular, for angle T12 (defined so that tu + Tqi + = ^) have 

^2^^2_^2 ^^^^ _ ^A(mf,mi,P) 



cos ri2 



sm ri2 



2mim2 



2mim2 

Here the "triangle" function A is defined as 

A(x, y, z) = 2xy + 2yz + 2zx — — y"^ — z^ = — A(x, y, z) 



(2.4) 



(2.5) 



where X{x,y,z) is the well-known Kallen function. The result (|2.2|) can be related to those 
presented in Ref. by a simple transformation of the occurring 2F1 functions. 
Using Kummer relations for the contiguous 2-^1 functions, we get 



1, e 



:i-2e)2Fiy I 



1 - 2£(1 - z) 2F1 



1, l + e 



(2.6) 



The resulting 2F1 function can be represented as 



I, l + e 
3 
2 



sin^^ 



sin ^ (cos ^) 1-2^ 



fe{e) , 



(2.7) 



with (see in 11, 



d0(cos0)2^ = 2-1-2^ J2 ^ [Ls,+i(^ - 29) - Ls,+i(7r)] 



(2. 



As a result, we reproduce the e-expansion of the two-point integral obtained in Ref. p4 

Til + e 



J(2)(4-2£;1,1) = ivr 



2-e 



mi^^+m2^^ 



2(1-2£)l £ 
[A(m2,m2,A;2)]'/'"" °° 



+ 



2 2 
— m2 



2e 



-m, 



-2s 



'-^ (2fV ^ ^ 
T.^—rT. [Ls,+i(vr) - Ls,+i(2r^J] . (2.9) 
1=0 J- 1=1 



(F)i- 



The expansion ( ^.9] ) is directly applicable in the region where A(mf,m2,/c^) > 0, i.e. when 
(mi — 7712)^ k"^ {ttli + 7712)^. To obtain results valid in the region A{m1,m2, k"^) < 
(i.e., X{m1,m2, k"^) > 0) the proper analytic continuation of the occurring Lsj(6') should be 
constructed. 

An important special case of Eqs. (|3)-(|]8D is e = ^ {z = \). It corresponds to the 
on-shell value fc^ = of the integral ( |2.9| ) with rrii = m2 = m. In this case, we obtain 
Lsj+i(^^j which correspond to the odd basis discussed in the introduction. Expansion of 

more general 2-^1 functions of argument 2 = | is discussed in Appendix B.2. 

We shall also need another representation for the two-point function (see, e.g., in ||3^ 



/2)(4_2£;1,1) = in^-'T{e) 



1 r^fi-^) 



F 1 r(2-2£) 



x'/^-%k~re 



2\e Aire 



m 



-2e 



2(1 -e) 

-2e 



+ 



2 1-e 



fc^ — + m\ — \J A(mf , mg, A;^) 



1, e 
2-e 



^2 



(2.10) 



where 



^1 



Am^k'^ 



Z2 



X{mi, ml, k"^) — m\ + m\ — k"^ 



Amfk'^ 



(2.11) 

with X{ml,ml, k"^) defined in Eq. (P3|). 

Employing Kummer's relations for contiguous functions, one can transform the 2-^1 func- 
tion from Eq. ( p.lOj ) into 



'1 + z-(1-z)SFi 



2(l-2e)^ 



f2.12) 



The resuhing 2-^1 function can be expressed in terms of a simple one-fold parametric integral, 

1 



1, l+e 



[l-z) 



Expanding the integrand in e, we get 

' 1, e 







[i-ztr-i 



(2.13) 



2-6 



1 - F \ 00 J ^ 

^^^^^^^l^l + z-{l-zy-'^-2il-zy-'^eY.e^^^^ 

(2.14) 



where Sa,b{z) is the Nielsen polylogarithm (see, e.g., in Ref. ||3^), whose definition ( |A.17| ) 
and some properties are collected in Appendix A.l. 



2.2 Analytic continuation 

To relate results in different regions, it is convenient to introduce the variables 

ln(— 2; — icrO) = ln{z) — icrvr. 



(2.15) 



where we put 9j = ^r^j, and the choice of the sign a = ±1 is related to the causal "+iO" 
prescription for the propagators. Whenever possible, we shall keep a undetermined, since 
one may need different signs in different situations. 

Let us note that, transforming from variable z to 1/z, we get the same 2-^1 function. 



2F1 



1,6 
2-6 



' r(2-L) '--''°(^") ■ f^-"' 
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In particular, we can replace each 2-^1 function in Eq. ( p.lO| ) by a linear combination of l.li.s. 
and r.h.s. of Eq. ( p.l6| ), with the sum of the corresponding coefficients equal to one (say, pi^2 
and (1 - pi,2)), 



J(2)(4_2£;l,l) 



. r(l +e) /m]" +m; 



2e I ^-2£ ^2 ^2 



ITT 



2{l-2e) \ e 



+ 



mi — mr, 
71^ 



+icr — 



Pi(-2i) ^ - (1 - Pi)(-2;j)' 



+2E 



i=l 



00 J 



j=i k=i 



00 J 



-(1 - p,)(-^.)^E^'' Y.(-2y~'Sk,-k+i{i/ 
j=i k=i 



(2.17) 



Comparing ( p.9|) and (|2.17| ) for pi = p2 = I, we arrive at the following analytic continu- 
ation of the functions involved in the e-expansion: 



ia [Ls,{n) - Ls,(^^)] = — ln^(-z) [l - (-l)^' 



+i-iyU - 1)! E E (-2)-' [Sk,-k-,i^) - i-irS,,.,.,{l/z)] . (2.18) 

p=0 ^ P- k=l 



where Zi and a are defined in ( |2.15|) . In particular. 



j=o J- 



Lsj+i(7r) - Lsj+i(6') 



1 

2e 



00 J 



-(-^)-^E^''E(-2r'^^,w(^) + (-^)^E^''E(-2r%w(v^) -(2.19) 

j=l k=l j=l k=l 

Since Lsi(^^) = —9, Ls2(7r) = we get 

ia[Lsi(7r) -Lsi(0)] = H-z), ia [Ls2(7r) -Ls2(0)] = -| [Li2(^)-Li2(l/2)] . (2.20) 

We note that for higher values of j the number of generalized polylogarithms involved in 
Eq. ( p.l8| ) can be reduced. For even j = 2/ (/ > 2), we have 

ia [Ls2z(7r) - Ls2i(6')] = 

(2/ - 1)! ln*'(-z) '"^"^ f, 
7^ E j — E cii,p+k[Sk+i,2i~k^p~i{z) - {-iySk+i,2i-k-p-i{'i/z)] 



p=o P- k=0 



k=0 



k\ 



Li2;-fc(^)-(-l)'Li2;_fc(l/^) 



(2.21) 
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whereas for odd j = 2/ + 1 we have 



where 



icr [Ls2«+i(7r) - Ls2i+i 
{2l)\ '"^ \if{-z) '"^"-^ ^ 



p=0 



k=0 



E M/, In (-^) Li2i-fc+i(z) -(-l)'=Li2i-fc+i(l/2;) 



42' /! ^t^^ 



^ '-^^ (p + 2g-l)! 



(2.22) 



(2.23) 



2.3 Massless limit 

As a simple example we can consider the limit when one of the masses vanishes, mi 
(m2 = m). Using hypergeometric representation (see, e.g., Eq. (10) of |4l]). 



/2)(4-2£;l,l) 



mi=0, m2=m 



r(i + ^) 

eil-e) 



2Fl 



1, e 
2-e 



4) (2-24) 



and Eq. ( p. 14 ), the result for an arbitrary term of the e-expansion can be obtained 

-2.r(i + £) 



j(2)(4-25;l,r 



mi=0, m2=m 



1 



X • 



1 -M r 
2ue 



'I 



u 



1 



'1 - 2e) 

j=i k=i 



u 



u 



(2.25) 



with u = k'^/m'^. 

Note that for this limit the terms up to order can be extracted from Eq. (A. 3) of 
Ref. Our expressions are in agreement with their results. We would like to mention, 

that each term of the expansion (|2.25| ) can be obtained from general results, (|2.9| ) or (|2.17|) . 
Expanding them with respect to ml and In(m^) (see details in [^), we reveal, that there is 
a limit mf = 0. This procedure allows to check the unambiguity of a: the limit 



must 

exist in each order of e. So, the finite part gives relation between ( 2.15 ) and ln(— m^/A;^), 
whereas linear one gives rise to correlation between sign of the analytical continuation of 
function \{ml,ml,k'^) and ln(— m^/fc^). Since, for Feynman propagator, ln(— m^/A;^) = 
ln(m^/A;^) + ivr, we have a = —1. 



3 Three-point examples 
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3.1 General remarks 

In this section we shall consider e-expansion of one-loop three-point integrals, shown in 
Fig. 1, 

J {n; ui, 1^2, 1^3) = / 77 2Tm7 — V2 21^^571 21^- (^•-'-) 

[{P2 - qr - ^1] [{pi + Q) - ^2] [Q - ^i] 




Figure 1: One- loop three-point function with masses rrii and momenta pi [pi + P2 + = 0) 



We shall mainly be interested in the case of positive integer powers of the propagators 
i/j. Using the integration-by-parts approach ^ (for details, see in |37]), all such integrals 
can be algebraically reduced to J^^\n; 1, 1, 1) and two-point integrals. Therefore, we shall 
concentrate on the case ui = 1/2 = I's = 1. 

To construct terms of the e-expansion of J^^\n; 1, 1, 1) with general masses and external 
momenta, the geometrical description seems to be rather instructive. The geometrical ap- 
proach to the three-point function is discussed in section V of (see also in [^). This 
function can be represented as an integral over a spherical (or hyperbolic) triangle, as shown 
in Fig. 6 of with a weight factor 1/ cos^~^^6' (see eqs. (3.38)-(3.39) of [^). This trian- 
gle 123 is split into three triangles 012, 023 and 031. Then, each of them is split into two 



rectangular triangles, according to Fig. 9 of [^. We consider the contribution of one of the 
six resulting triangles, namely the left rectangular triangle in Fig. 9. Its angle at the vertex 
is denoted as |v'^2 5 whereas the height dropped from the vertex is denoted 7712. 
The remaining angular integration is (see eq. (5.16) of 



1 

2e 



dip 



1 + 



tan^?7i2\ 

COS^ (f J 



IE 



d(f In-' 



+1 



14 



tan"' r]i2 



cos^ ip 



(3.2) 
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First of all, we note that the l.h.s. of Eq. ( p.2| ) yields a representation valid for an arbitrary 
e (i.e., in any dimension). To get the result for the general three-point function, we need to 
consider a sum of six such integrals. The resulting representation is closely related to the 
representation in terms of hypergeometric functions of two arguments (see also in |^| 
for some special cases). 

In the limit £ ^ we get a combination of CI2 functions^, eq. (5.17) of p2|. Collecting 
the results for all six triangles, we get the result for the three-point function with arbitrary 
masses and external momenta, corresponding (at e = 0) to the analytic continuation of the 



well-known formula presented in The higher terms of the e-expansion correspond to 

the angular integrals on the r.h.s. of Eq. (|3.2|). We note that the e-term of the three-point 



function with general masses has been calculated in |^ in terms of Lis. 

An important special case is when all internal masses are equal to zero, whereas the 
external momenta are off shell. This case was considered in detail in Refs. |24, 25]. The 
following result was obtained: 



j(3)(4-2£; 1,1,1) 



mi=m2=m3=0 



2^2-. .1+2. ni+e)V\l-e) [^{pIpIpI 



2 J2 ^2M-l/2+£ 



r(l-2e) 



^ (-26)^ 



Lsj+2(7r) - [Lsj+2(vr)-Lsj-+2(20i)] 



i=l 



where the angles 0j {i = 1,2,3) are defined via 

, pI+pI-pI , pI + pI-pI , pI + pI-pI 

cos 01 = p== , COS 02 = p== , COS 03 — 



,(3.3) 



^\^p¥i 



2\^Pipi 



'i^pipi 



(3.4) 



SO that 01 + 02 + 03 = TT. Therefore, the angles 0, can be understood as the angles of a 
triangle whose sides are ypf , -y/pf and whereas its area is \\l ^{piipiiPi)- Introducing 
Oi = 20j, we can use the same procedure of analytic continuation as discussed in Section 2.2 
(with cr = 1), in terms of the Nielsen polylogarithms Sa^biz). 

Note that the integral (|3.3|) is closely related to the two-loop vacuum integral with dif- 



ferent masses, due to the magic connection The result for this integral, which is called 
I{n; z/i, i>2, t's), is also presented in Refs. p^ , pS] , in an arbitrary order in e. For the analytic 
continuation, we can also apply the procedure of Section 2.2 (with a = —1). Such integrals 
will also occur in Section 4. 



3.2 On-shell triangle with two different masses 

Consider a three-point integral ( p.l| ) with one massless propagator (7713 = 0) and two adjacent 
legs on shell (pg = ""^i. Pi = ^V)- Such integrals are important, e.g., for studying corrections 
to the the muon decay in the Fermi model. Using Feynman parametric representation, one 

^°For e = i (n = 3) we reproduce the well-known result of in terms of elementary functions (for 
further details, see Section VA of ||2^). 
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can show that in the z/i = 1/2 = ^3 = 1 case this integral can be reduced to a two-point integral 
with the masses of internal particles mi and 1712 and the external momentum k = p^, 



j(3)(4-2£; 1,1,1; 



71 



rn^=0, p\=rn\, p\=rn'^ 2£ 



j(2)(2-2£;l,l) . 



(3.5) 



Note that all powers of propagators are equal to one, whereas the space-time dimension is 
(2 — 2e) in the two-point integral. 

Using Eq. (6) of we can represent it in terms of the (4 — 2£) -dimensional integrals. 



J(^)(2-2£;1, 1) 



TT 



(4 - 2e; 2, 1) + (4 - 2e; 1, 2] 



(3.6) 



Then, using the integration- by-parts technique p[ (see also Eq. (A. 17) of [0), we can express 
the integrals with second power of one of the propagators in terms of J'^^)(4 — 25; 1, 1) and 
tadpole integrals. As a result, we get 



j(2)(2-2e;l, 1) 



2(l-2£)A;2j(2)(4-2£;l,l) 



-i TT^"^ r(e) Uk"^ + ml- ml)m^'^^ + (A;^ -ml + mDm'^'^'' 



(3.7) 



with A(mi,m|, k"^) defined in Eq. ( p.5|) . 

In particular, the e-expansion of this integral is 



j(^)(2-2.; 1, 1) = -i vr^- ^ . + gf^,^^ E ^ E i^^M^) ' Ls,+i(2r-J] . (3.8) 



[A(mf , 777,2, ^^)] 



j=0 



This gives us the e-expansion of the on-shell triangle diagram considered. Extracting the 
infrared (on-shell) singularity, it can be presented as 



J(=')(4- 25; 1,1,1) 



ivr^-^ Til + e) 



m3=0, p'^=m\, p\=m'^ 

00 2 



{k 



2\e 



[A(mf , 7772, k'^)] 



1/2+e 



x^^-E 



{2e) 



E[Lsj+2(7r) -Lsj+2(2%)] 



(3.9) 



where ri2 = vr — Tq^ — rQ2. 

Since the functions are of the same type as in the two-point example, we can use Eq. ( |2.18| ) 
for the analytic continuation. 



3.3 A more complicated example 



Let us consider triangle integral (|3.1| ) with 7771 = 7772 = 7773 = 777, = ^2 = 0, with an 
arbitrary (off-shell) value of p\. Such diagrams occur, for example, in Higgs decay into two 
photons or two gluons via a massive quark loop. Following the notation of Ref. 



we 
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shall denote this integral (with unit powers of propagators) as J3(l,l,l;m). According to 



Eq. (40) of Ref . , the result in an arbitrary space-time dimension n = 4 — 2e is 



Jfifl, 1, 1; m 



5 -^5 -^1 "V \pf=p^=0 



Am? 



(3.10) 



For the occurring 3F2 function, various one-fold integral representations can be con- 
structed (see, e.g., in [^). Below we hst some of them: 



Vil-e) 



2zr(i + £)r(i-2e) 7 t^i^ 



ln(l - tz) 



4(1 -t) 



— / dr 

ze J 


1 



1 - 



cosh 



- 1 



^/^{[i-4..(i-.)r-i}. 



(3.11) 
(3.12) 
(3.13) 



In fact, to get representations of the terms of the ^-expansion in terms of Clausen and 
log-sine functions, the following two-fold angular integral representation appears to be rather 
convenient: 



sin^ e ,K 



1,1, l + £ 

3-^2 I 3 2 



It is easy to see that 



sin^ej = 2 j d0(cos0)-^"y d0(cos. 
^0 



\2e 



Y.e'h,{e) (3.14) 



hi{9) = 2CI3 (tt - 2^) - 2013 (tt) - 2^Cl2 (vr - 2^) 



(3.15) 
(3.16) 



Moreover, in terms of the function whose definition and the e-expansion are given in 
Eq. ( |2.8|) , we can represent the integral in Eq. (|3.14|) as 



2 J d0(cos0)-^"y d0(cos. 




\2e 



2 J fM d/_,(0) = 2fMf-eiO) - 2 j d/.(0) 




fe{e)f-e{e) + J [fM d/_,(0) - f.M d/.(0)] • (s.ir) 




The last representation is nothing but splitting the function ( |3.14| ) into the even and odd 
parts. The first term, f^{d)f-e{d)i is symmetric with respect to e ^ — e and therefore 
contains only even powers of e. The second (integral) term is antisymmetric and contains 
only odd powers. 
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Thus we have obtained results for all hj{9) with even values of j (j = 21). One needs just 
to pick up the e^' term of the expansion of /e(6')/_e(6'), with given explicitly in Eq. (^78|) . 
For example, 

h2{9) = -[Ls2(7r-2^)]' + 2e[Ls3(7r-2^) -Ls3(7r)] , (3.18) 
h^{e) = [Ls3(7r-2e)-Ls3(7r)]^-|Ls2(7r-2^)[Ls4(7r-2^)-Ls4(7r)] 

+|0[Ls5(7r-20) -Ls5(7r)] . (3.19) 

However, the calculation of the odd terms of the expansion (j = 2/ + 1) is less trivial, 
starting from j = 3 (the result for hi{9) is given in Eq. (p.l5|) ). It appears that hs{9) (and 
the higher odd functions) cannot be expressed in terms of the Lsj functions and their simple 
generalizations. 

Of course, for a given value of 9 there is no problem to calculate hj{9) with very high 
precision, using integral representation ( |3.13| ). As an illustration, let us consider a special 
value of pI, pi = rn? {z = \, 9 = |). Then we have 



(f ) 



1 



2(J + 1) 



— In^+Vl-x + x^) 

X 



(3.20) 



In particular, for h^{^) we obtain a number which, with the help of the PSLQ program p9 
can be identified as 



^3(1 ) - + ^(2(3 + iiCs - ^7rLs4(|) + ^7rC2Ls2(f) 



(3.21) 



where 



X5 = S?;- 2(1) 



°° {n\? 1 



E 



°° 'n\Y 1 



[n) + 7]= 



„tl (2n)! ' {2n)\ 
0.0678269619272092908692628002256569482253389 . . . (3.22) 



is a special case of binomial sum |17, 33]. Therefore, the obtained e-expansion of the 3F2 



function at 2; = | is 



3-^2 



1,1, l+£ 
^ 2 

2' ^ 



C2 + £^fC3-f7rLs2(f 



e^|lOC4 + |7rLs3(f)-f [Ls2(f) 



+^'{^5 + fM^ + |^C5 - l^Ls4(f ) + ^VrC2LS2 
Ls3(f)]'-fLs2(f)Ls4(f)+|vrLs5(f 



+47rC2Ls3(f ) + f 7rC3Ls2(f) + ^^C,^+0{e 



(3.23) 



In Appendix B we present several terms of the ^-expansion of the 3F2 function with more 
general parameters. 
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For another special case 



3m^ 



l,e = f ), we find 



7rC2LS2 



\tx\jSa 



+ tC2C3 + ICs 



(3.24) 



Note that here we do not need that extra basis element xs- 

Let us discuss the appearance of the constant X5 ( p.22| ) and its relation to transcendental 
numbers examined in |^ . A detailed analysis of hyper geometric functions of argument z = ^ 
(see Appendix B) and Feynman diagrams (see Eqs. ( p.lO| ), ( ^3|) and ( [4.10| )) allows us to state 
that at the level 5 there is a new (with respect to the constants defined in |2^) independent 
irrationality related to the binomial sum ( p.22|) . The set of sums or their linear combinations 
involving this new irrationality (for a more detailed discussion, see Appendix B.l) consists 
of six elements: S^lJ.g, S^lJ.g, Si'I.3 and three linear combinations given in Appendix B, 
Eq. ( p.7| ). For example. 



X5 + 2S?;- 3(1) = |7rLs4(f ) - f C2C3 - f Cs. 



Of course, instead of ( ^.221 ) one could choose any other element of this set. The linear 
independence of xs (|3.22| ) and other weight-5 elements of the odd basis was established by 
PSLQ-analysis of all these constants calculated with 800-decimal accuracyj^. 

Introduction of a new constant ( p. 221 ) means that the ansatz for the odd basis of weight 5 
elaborated in Ref. |2^ (see Table 1 of [28||) should be re-analyzed. We note that one of the 
elements, :^Cl5 is proportional to :^C5- Although this element is linearly independent, 
its appearance contradicts one of the mnemonic principles of constructing the ansatz: any 
non-factorizable element may appear only once, with or without the normalization factor 
Since we already have Cs (without -^), this element can be omitted. Including, instead 
of it, the new element xs? we are able to save the rule Nj = 2\ whereas the simplicity of the 
ansatz gets spoiled, since, as far as we know, the new element is not expressible in terms of 
the Clausen or log-sine functions of | or 

Let us also consider cases when hj{9) can be expressed in terms of the even basis. The 
threshold case, p| = 4m^ (z = 1, 6' = |), is trivial: 



3-^2 




dA 



1-A 



-2e 



A2 



1 r 

2e 



-i 00 



i=o 



Therefore, we get 



(f ) = 1(2^'^^ - 1)C 



•i+2- 



ljO+2. 

(3.25) 
(3.26) 

In particular, = -yCs- 

Another case of interest is = t- Here, the situation is rather similar to the case 6 = ^. 

Again, the calculation of /13 (^jj happens to be non-trivial. With the help of the PSLQ 

^^The algorithm for high precision calculation of hsj{9) and Lsj^''(0) functions is discussed in Ap- 
pendix A of Ref. |3^ . The corresponding programs can be found on Oleg Veretin's home page at 



http : //www . if h . de/T veretin/MPFUN/mpf un . html 
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program, we obtain 

hs (f ) = 1^X5 - fvrCU (f ) + ^nCCh (f ) + - §-,(2(3 , (3.27) 

where 

00 |'„n2 on 00 („]\2 on 

~ 0.4951660770553611208168992757694958708378.... (3.28) 

One can see that this multiple binomial sum is nothing but the k = 2 version of the odd 
element Xs given in Eq. ( |3.22D . Therefore, this element Xs should be added to the set of 
weight-5 elements of the even basis. A further discussion of the even-basis elements is given 
in Appendix B.l. 



4 Two- and three-loop examples 
4.1 General remarks 

Below we present the e-expansion of several master integrals arising in FORM packages 
for calculation of two-loop on-shell self-energy diagrams and three-loop vacuum integrals 



31, Our method of calculation is a combination of the Mellin-Barnes technique 



and the PSLQ analysis [^: we present results for the master integrals in terms of 



hypergeometric functions, whose expansions, given in Appendix B, are obtained by means of 
PSLQ. In some cases we give the result for arbitrary mass and/or momentuml^. Sometimes, 
the integrals with higher powers of some propagators happen to be simpler from the e- 
expansion point of view (for instance, the integrals having an ultraviolet-divergent subgraph, 
like VlOOl or JOll in Fig. H). In such cases we derive an expansion of these "simple" 



integrals and then reduce them to the master integrals by means of packages [0] or ||3T| . 

The integrals under consideration and notations (indices and mass distributions) are 
shown in Fig. |^. We are working in Minkowski space-time with dimension n = 4 — 2e. 
Moreover, each loop is multiplied by the normalization factor [i7r"''^r(l -|-e)]~^. For example, 
the two-loop vacuum integrals with different masses, which were denoted in Refs. |1^, ^ 
as I{n; Ui, 1/2, 1^3; rrii, 1712, fris), are in the case of equal masses normalized as 

I{n;ui,U2,J^3;m,m,m)\^^^ = -7r'T^(l + £:) VLlll(i^i, z/2, z/3) (4.1) 

where we adopt the notation VLlll from Ref. pHfQ (and also put m = 1). 

We can mention an example of physical calculations |^3|, where the finite parts of E3, 



D3 (which was found in p8|, BOl) and D4 (found in Ref. |T^) have been used 



"'^^A review of modern computer packages for calculations in high-energy physics is given in p9| . Some 
useful details concerning recurrence relations for multiloop integrals can be found in |5C| | . 

""^•^We would like to mention that in Ref. pl[ divergent parts of some three-loop vacuum integrals with two 
different masses have been calculated, for a special case when the ratio of the masses is \/3. 

^''Two-loop bubble integrals are also implemented in the package js^ . 
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Vl 



D4 (Vi,Vj,V3,v^,V5,Vj) E3 {cTi,a,,Pi,Pj,v) D3 (0,a,0,Vi .v^.Vj) 

Figure 2: Two- and three-loop diagrams considered in the paper. Bold and thin lines 
correspond to massive and massless propagators, respectively. 



4.2 FlOlOl 



This integral is a good illustration of the application of general expressions given in Ap- 
pendix B. The off-shell result for this integral in arbitrary dimension was presented in p4| 
(where it was called /s, see Eq. (22) of For unit powers of propagators, the result reads 



2^^)Flolol(^»^m) 



l,l + e,l + £,l + 2£ 
1 + e,2 + e,2- e 
2 \ 




4^2 



1,1, l + £ 



P 



Am? 



(4.2) 



The finite part of -Fioioi is given in |5^, ^ , whereas the term linear in e has been considered 



m 



33 1 . An algorithm for the small momentum expansion of such diagrams (involving a "zero- 



threshold" ) has been constructed in Ref. where also explicit results for a few terms have 
been given, for general values of the masses of massive particles, in the limit £ — 0. All 
these calculations agree with the result ( |4.2| ). 

For p2 > this diagram has imaginary part, due to the factor (— m^/p^)^, which is related 
to the massless two-particle cut. If < 4m^, this is the only contribution to the imaginary 
part, whereas for > Arn^ we get additional contributions from the cuts involving two 
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massive particles. Therefore, for < < Aiv? we get 



2+ie, 



1 -2e)ImFioioi(r,m) = - 



r(l-e) 



.2\ ^ 



2T{l + e)T{l-2e) \p^ 



^1 



1,1, 1 + e 
3 
2 



3- 2 

> 5 ^ 



4m2 



(4.3) 

with the same 3F2 function as in the three-point example in section 3.3. This is expected, 
since, cutting the considered two-loop diagram across the two massless lines, we get nothing 
but that triangle diagram. Therefore, all results for the e-expansion {hj{6) functions, etc.) 
presented in section 3.4 are directly applicable to the imaginary part of Fioioi- 

In particular, in the on-shell limit p"^ = m? (we also put m = 1) we follow the notation 
of Ref. 1^ and denote 



FlOlOl = Fi 



10101 



p^=m?; m=l 



(4.4) 



where it is understood that all powers of all propagators are equal to one. This integral is 
used as one of the master integrals in the 0NSHELL2 package |^ . 

Using expressions given in Appendix B, we can now present the ^-expansion of the real 
part up to the e"^ term. 



(1 -2£)ReF10101 



-4C3 + 27rLs2 f 



'|67rLs4(f 



Ls2(f)]^-7vrLs3(f)-fO 



ix5 



^C2C3 - IC5 + ivrC2Ls2(f ) - ivrLs4(f ) + 0(e^)(4.5) 



The finite and linear in e parts coincide with those given in |^ and |Q, respectively^^. 
According to Eq. (O), the imaginary part can be extracted from Eq. (|3.23| ), 



;i - 2e)Im FlOlOl 



TT 



2r(i + £)r(i-2e) 



[r.h.s. of Eq. (3.23)] , 



(4.6) 



which yields its e-expansion up to e'^ 



4.3 D, 



The integral 04(1, 1, 1, 1, 1, 1) is another master integral used in ||31|, The finite part of 
this integral was given in . Here we are going to obtain result for the e-term. It is easy to 
see that the integral 04(1, 1, 1, 1, 1, i/) can be obtained by integrating the off-shell two- loop 
diagram Fioioi (see ||5^ Eg. (|4.2| )) with a massive propagator raised to a power z/. Using 
the following property (in our case = Am?) 



1 



ITT 



n/2 



p 



(p2 _ ]Vf2)i 



( ai,- 


■ ,ap 


P'\ 


^ &i,- 







""^^ We take into account the Erratum to J27 
there is another misprint in Table I of Ref. 
«! and a2 should be changed: ai = — 1, 02 = ; 
in (H). 



where Eq. (1) and Table 1 of |^ were corrected. In fact, 
for the master integral Flllll the signs of the coefficients 
The corrected result is presented in p. 541 of |^ (see also 
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X p+iFq+i 



r(z.-f) 



r(t-z.) 



Q+1 



r(t) 



Q j^ r(&,)r(z. + a.-f) 



1^. 



ai - |, . . . , z/ + ap 



2 ' 



V 



(4.7) 



where {aj, 6j} ^ 0, —1, —2, . . ., we arrive at a resuh which contains (for an arbitrary v) three 
4F3 functions and three 3F2 functions. For z/ = 1, two of the three 4F3 functions reduce to 
the 3F2 functions of the type considered in Appendix B. For some 3F2 functions we use the 
the relation 



ai, ■ ■ ■ , ap 



z- 



ai . . . ap 



1, 1 + ai, 
2,1 + 61, 



, 1 + ap 



(4. 



Then the master integral becomes 
(1-£)(1-2£)D4(1, 1,1, 1,1,1): 

1 



' e{l-e){l + e){l + 2e) 



4-^3 



l,l + £,l + £,l + 2e 
3 
2 



f+£,2 + e,2-e 



r(l-e)r2(l + 2£)r(l + 3e) 
'e{l + 2e){l + 4£)r2(l + e)r(l + Ae) 



I l,l + 2e,l + 2e,l + 3e 
^ M I + 2e,2,2 + 2e 



1 



2£2(l + £^ 



3-^2 



l,l + £, 1 + e 
|,2 + £ 



r(i-£)r(i + 2£) 



3-^2 



1,1, l + £ 

^ 2 

2' ^ 



1 



2^2(1 + 2e)2 



3-^2 



1, 1 + 2e, 1 + 2e 
I + e, 2 + 2£ 



4e2r(l +£) 

r(i-£)r2(i + 2£)r(i + 3£) 



4e4r2(l + £)r(l + 4e) 



r(i -£)r2(i + 2£)r(i + 3£) 

^4^2(1 + 2e){l + 4£)r2(l + e)r(l + Ae) 



3-^2 



1, 1 + 2£, 1 + 3£ 
I + 2£, 2 + 2£ 



1 

4? 



(4.9) 



Now all pFq functions belong to the types considered in Appendix B. Using those results, 
we obtain the following result for the lowest terms of the ^-expansion: 



;i-£)(i- 25)04(1,1,1,1, 1,1) 



2C3 



+^ {f X5 - 1^C2C3 + If C5 + 1 vrC2Ls2(f ) - l^Ls4(f ) + 7vrLs4(f )} + 0{e'). (4.10) 



The finite term is in full agreement with the result of |13[ , whereas the result for the e-term 
is new. As we see, the new constant X5, Eq- ( |3.22| ), also appears in the e-part of this integral. 



4.4 VlOOl 



The off-shell integral Viooi with two different masses and unit powers of propagators was con- 
sidered in The result, Eq. (46) of |^ is presented as a sum of F2 and F4 hypergeometric 
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functions of two variables 
Again, we define 

V1001(a,cri,a2,/?) = Viooi(a,o"i,a2,/3;p^,m,m) (4.11) 
The finite part of on-sliell master integral has been calculated in , while the linear in e 



term is presented in |2^. The on-shell diagram V1001(l, 1, 1, 1) belongs to the set of master 
integrals used in the 0NSHELL2 package Employing the Mellin-Barnes technique |^ 



one can find the following result for an arbitrary set of the indices (remember that we put 
m = 1): 



V1001(a,ai,(T2,/3) 



r(f-ai)r(f-a2)r(cri + a2- 



r(a + ;3 + (Ti + (T2 - n)V{2n - 2/3 - a - 2ai 



X 



2^2) 



r (^^n — a — (3 — Gi — 02 

^ / (3,a + (3 + ai + a2 - n,l + a + (3 + ai + a2 - \' 
^ ^ ^ V I + ^ + (3 + ai + a2 - n,\{a + I) + 13 + ai + (72 - n 

r (f ) r (f + /? + di + (T2 - n) r (n - - (72 
2r(/3)r 



X 



3-^2 



- 1 + - 

2 ' ^ 2 



2,n - (Tl 

il + n- /9-(Ti- 



- (T2 - ; 

0^2-9 



r (^) r (^^^ + /5 + (Ti + (T2 - r - ai - a2 - ^ 



0-1 



X 3^2 



2T{I3)T 

\{a + l),\{?> + a — n),n — Gi — 02 + \{l — a) 



f + n 



(4.12) 



2 ' 

^ 14-/^-^1-^2 

Let us consider the case di = ct2 = /3 = 1. For a = 1 (the master integral), and also for 
a = 2, the result ( [4.12 ) can be expressed in terms of the 2F1 functions. For a = 1, they 
can be transformed using the relation ( p.6|) and another combination of Kummer relations 
for contiguous functions, 

2a -3e)2 ^ 12^(1 - 3e);z(l 




25) 



(l-25)(l + 4e) 




(4.13) 



The integral with a = 2 can be reduced by the 0NSHELL2 package to this master 
integral, plus some combination of F-functions: 

F(l-£)F(l + 2£) 



V1001(2, 1,1,1) 



-i(l-2£)V1001(l, 1,1,1) + 



3^2(1 - 2£)F(1 + £) 



(1 - Ae)V^{l - £)F(1 - 4£)F(1 + 2e) 
'6^2(1 - 3e)(l - 2£)F(1 - 2£)F(1 - 3£)F(1 + e) 



(4.14) 



^^We note that the occurring function can be reduced to 2-F1 function of the type (2.7), whereas the 
F2 function can be transformed into Fi function of two variables (see Eq. (A. 64) in |60| ). 
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This connection provides a non-trivial check on the results for the e-expansion of the occur- 
ring integrals. The result for the master integral can be presented as 

3T{1- e)T{l + 2e) 



;i-2£)^V1001(l, 1,1,1) 



8eT{l + e) 




l,l + £ 



8e{l + 4£)r(l - 2£)r(l - 3£)r(l + e 
(1 - e)T\l - e)T{l - Ae)T{l + 2e) 



4e2(i - 3e)r(l - 2e)T{l - 3£)r(l + e 
£)r(l-4£)r(l + 4£) 




1,1 + 3£ 
+ 2e 



6)T{l + 26) 



7r35-T3(l 



2eT^{l - 2e)T{l + 2e)T{l + e) 



(4.15) 



Here we have two 2-^1 hypergeometric functions of the argument ^. The ^-expansion of the 
first function can be extracted from |Q (see also Eqs. ( |2.7| )-( P^ ) of this paper). 




E 



{2ey 



7 32-^—0 
whereas the second function can be represented as 



Ls 



■j+i 



(f ) - L^.- 



(4.16) 



2F1 



l,l + 3£ 

l + 2e 



2^+6^(1 + 4e) 



7r/6 



3t 



+e 



d0(sin( 



\4e. 



COS( 



\2e 



(4.17) 



When we expand in e we get integrals of the products of various powers of ln(sin0) and 
ln(cos0). An integral representation of this function for an arbitrary argument is given 
in Eq. ( p.l3|) . Its e-expansion can be written in terms of Lsc function (for details, see 
Appendix A. 2). The expansion up to looks like 



3t 



+£ 




-TX 



(2£)|Ls2(f) + {2ef f^TTCa - iLsgff 



1152 



Cs - iLs4(f) - |LS4(^ 

C4 + !C2Ls3(f ) - iLs^d) - ^Ls,(f ) - |vrLsi^)(f ) + fLs?) (f 



4(1 + As) 
+ (2£)3 

+ (2^)' 

+0{e^) . (4.18; 

These expressions allow us to fix an error in the e part, presented in Eq. (2) of f^^'' 
The correct expression is: 



V1001(l, 1,1,1) 



2e 



1 1 
+ - 




+ IC2 



TT 



7- 



V3 V3 



+ ^ln3 



^"^ In Eq. (2) and Table 2 of |2^], results for two integrals, Vllll(l, 1, 1, 1) (with all massive lines) 
and V1001(l, 1, 1, 1), are combined. The result for Vllll(l, 1, 1, 1) is correct. In order to correct the 
result for the e-term of V1001(l, 1, 1, 1), the following changes should be done: (i) the term +952-^ (with 

5'2 = ^Ls2(f )) should read +|6452^, where 64 = for V1001(l, 1, 1, 1), according to Table 2; (ii) the 

coefficient 63 of the Cs term (given in Table 1) should be rather than ^. The result for V1001(l, 1, 1, 1) 
on p. 546 of 1^ should also be corrected, according to Eq. ( 4.19 ) 
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+e|f + 8G + fCs - 12-^ - 28—^ + 7—^ ln3 

1 TT . o 21 TT . 2lLS3^2^ 



Our expressions (|4.15|) -( fi.l8|) provide two new terms, £^ and e^, of the ^-expansion of 
V1001(l, 1,1,1). 



4.5 E3 

The three-loop vacuum integral E3 with all indices equal to one is one of the master integrals 



used in Avdeev's package and in MATAD The result for general values of the indices 
is 

r(f -ai)r(f -cT2)r(ai + a2-f) 



E3(cri,Cr2,pi,p2,Z^) 



r(ai)r(a2)r(f)r3(3-|; 



X 



[ r(pi + p2-f)r(z/ + 0^1 + ^2- 
1 r(z/)r(pi + p2) 

r(z/+pi+(Ti + a2-n)r(z/+p2 + cri + cr2-n)r (z/+pi+p2 + cri + cr2- ^) r(n-z/-cri-(T2) 



^[ Pi, P2, Pi + P2 - f ,^ - 0-1 - 0-2 

^ H |(Pl+P2), 5(Pl+P2 + l),'^-'^-0"l-Cr2+l 



r(pi)r(p2)r(2z/ + pi + P2 + 2(Ti + 2(T2 - 2n)r(n - cTi - (72) 



3n 



/ I^, + Pi + (Ji + (72 - + P2 + O"! + 0-2 - ^, + Pi + P2 + O"! + 0-2 - ^ 

^ y Z/ + cri + (T2 + |(pi+p2)-n, Z/ + (7i+0-2 + ^(Pl+P2 + l)-"', '^ + Cri + Cr2-^ + l 



(4.20) 



For the purpose of expanding in e, the integral £3(1, 1, 1, 2, 1) is the simplest one. 



£3(1,1,1,2,1) 



1 



r(i-5)r(i + 2£) [ /1,1 + e 

Ae^{l-e){l-2e) T{l + e) \' | 



2r(i + 2£)r(i + 3£) / i,3£ 



(4.21) 



The £-expansion of the first 2-^1 function is given in Eq. ( ^.16| ), whereas the second one can 
be reduced to (|4.18| ) via Kummer relation 



2F, 



l,3e 
\ + 2e 



l + 4e 



i + 2e 



(4.22) 



which is a particular case of ( [4.8|) . 

Then, using recurrence relations PT|, E3(l, 1, 1, 2, 1) can be related to to the master 
integral E3(l, 1, 1, 1, 1) as 

r(i-£)r(i + 2£) 



(1-2£)E3(1,1,1,1,1) 



-3E3(1,1,1,2,1) - ■ 

(1 - 4£)r(i - £)r^(i + 2£)r(i + 3^) 

"3^3(1 -e){l- 2e){l - 3e)T{l + Ae)T^{l + e) 



(4.23) 
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In this way we get 
(1-2£)2E3(1,1,1,1,1) 



3r(i-e)r(i + 2e) 



3 
2 



4£2(i -£)r(l + e 

3r(i -£)r2(i + 2£)r(i + 3£) 



l,l + 3£ 
1 + 25 



4£2(l - + 45)r(l + 4£)r2(l + £ 
r(l - e)T{l + 2e) r(l - e)T\l + 2e)T{l + 3e) 



2e3(i _ e)r(l + e) 6e3(i - 3£)r(l + 4e)r2(l + e) 



i) 

(4.24) 



with the same 2F1 functions as in the case of VlOOl. The e-expansion of these 2F1 functions 
is given in ( ^.161 ) and ( [4.1(^ ). Thus, for this integral we have also produced two new {e^ and 
e^) terms of the ^-expansion, reaching the level of 6- loop calculations (see in [Q). We also 
find an error in the e-part of the result given in Eq. (10) of which should read (remember 
that we put m = 1 and Ls2(|) = ^^82) 



E3(l,l,l,l,l)--|,-^ + iri4 



Ls2(f) 



^ V3 



+ 



+J-^ + 102^''^=^ 



Ls^d) _^Ls2(f 



^^ln3-5C2-K3 + t^ + 9- r- 



+45 



Ls3(f 



9 



Lssff 



30- 



LS2 f 



ln3 + 3 



In 3 + 



.Mi) 



Ls2(f) 

LS4f^ 



ln2 3 



I 25 ^ A 5 ^ ^ q 13^ 94 ^ ^ 

+ y^S2 - 3— C,2lnd - —(,3 - — — 



' 73 

TT 

9 



17C2 



LS2 f 



(4.25) 



4.6 D3 and JOll 

The off-shell result for the sunset-type integral Jon with arbitrary powers of propagators has 



been obtained in by using the Mellin-Barnes technique 

r(z/i+z/2+a-n)r (f-o-) r (z/2+a-f ) r (z/i+^-f 



Jqii{ct,i'i,-U2]P ,m) = (m ) 



X4F3 



2\n—a-—ui—U2 



r(z/i)r(z/2)r f r(z/i + z/2 + 2a - n)T^ 3 - f 



a, ui + iy2 + cr - n, iy2 + a 



2 ' 



1^1 + (T 



^ a + \{ui + U2 - n) , a + \{ui + U2 + I - n) Am? 
The three-loop vacuum integral D3 can be obtained by integrating ( ^4.26| ) over p 

1 r d"p 



D3(0,cr,0, Z/i,Z/2,Z/3) 



iW2r(3-f] J (p2-m2l'^3 



m=l 



.(4.26) 



(4.27) 
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In such a way we find 
03(0, a, 0, z/i,z/2,i/3) 



X • 



r(z/0r(z/2)r(z/3)r (f ) r3 (3 - f ) 
r (z/i + cT - f ) r (z/2 + a - f ) r(z/i + + ^ - n)r (z/3 - f ) 



X aF:- 



r(z/i + z/2 + 2cr - n) 
cr, z/i + (T - ^, z/2 + (J - ^, z/i + z/2 + (T - n 



^ ^ V t - ^^3 + 1, + K'^i + ^^2 - (X + |(z/i + z/2 - n + 1) 
r (f - z/3) r (z/3 + a - f ) r(z/i + z/3 + a-n)r(z/2 + z/3 + a-n)r(z/i + z/2 + z/3 + a- f ) 



r((7)r(z/i + z/2 + 2z/3 + 2(7 - 2n) 

z/3 + a - ^, z/i + z/3 + a - ra, z/2 + z/3 + a - n, z/i + z/2 + z/3 + a - 



3n 



X 4i^3 



Z/3 - f + 1, Z/3 + i(z/i + Z/2) + (T - n, Z/3 + i(z/i + z/2 + 1) + 



(T — n 



-. (4.28) 



The integral D3(0, 1, 0, 1, 1, 1) is also one of the master integrals used in the packages 



The result including the e-term is available in |^ and |3^. Here we are going to provide 
some further terms of the ^-expansion. 



Again, when using the general result ( [4.28|) it is simpler to consider an integral with 
shifted indices, namely D3(0, 1, 0, 2, 2, 2), which is related to the master one through a simple 



relation 



9D3(0, 1, 0, 2, 2, 2) = -^TT^ + 2(1 - 3£)(1 - 2£)(2 - 3^)03(0, 1, 0, 1, 1, 1) . (4.29) 
Considering Eq. ( |4.28| ) at cr = 1, z/i = z/2 = z/3 = 2, we get 



:i-£)D3(0,l,0,2,2,2) 



3-^2 



l,l + e,l + 2e 
^ +e,l-6 



6(1 + 2e) ^2 
1 T{1- £)T^{l + 2e)T{l + 3e) 



9 

1 + 2£, 1 + 3£ 
l + 2e 



e{l + Ae) r2(l+£)r(l + 4£) 
Moreover, the occurring 2F1 function can be reduced to a product of F-functions^^ 



- .(4.30) 



2F1 



l + 2e,l + 3e 
| + 2e 



27r 



r(2 + Ae) 



A) 3l+3^r(l + 2e)r2(l + e) 



(4.31) 



see (|B.14|) . The 3F2 function in Eq. (^4.301) belongs to one of the types considered in Ap- 
pendix B. Its e-expansion is given by 



/ l,l + £,l + 2£ 

3^2 |+e,l-5 



2s 



35 



i+3£ 3 



+ 4£Ls2(f) 



""^^Eq. ( 4.31 ) can be obtained via a simple transformation of Eq. (31) on p. 495 of |Q. A standard way to 
obtain formulae of such type is to use Bayley's cubic transformations of hypergeometric functions. 



25 



+6" 



fvrC2 + 18Ls3(f) 
216C2Ls3(f) 



112^^3 + |Ls4(f)-36Ls4(f 



19^^, _ If Ls5(f ) + 54Ls5(f ) - 1087rLsl^)(f ) + SlLsf (f ) 



Finally, we obtain for the master integral: 
2e{l - e){l - 2e){l - 3£)(2 - 3£)D3(0, 1, 0, 1, 1, 1 



(4.32) 



9 



3-^2 



l,l + £,l + 2e 

§+£,!-£ 



(1 + 25) 

(4-15^) 67r r(l + 2£)r(l -£)r(l + 3£) 



32 



h3e 



14(1+ 



(4.33) 



Considering result (|4.28| ) with a different set of indices, cr = z/i = 1 and 1/2 = = 2, we 



get 



03(0,1,0,1,2,2) 



1 



'l-e)\l + 2e)e 



3-^2 



l,l + e,l + 2e 
I + £, 2 - e 



1 r2(i + 2£)r(i -£)r(i + 3e) 



V{l + Ae)V^{l + e) 



2e,?>e 
\ + 2e 



. (4.34) 



On the other hand, using the recurrence package |3T| we obtain a much simpler result for 
the same integral, 

D3(0, 1, 0, 1, 2, 2) = ■ (4.35) 

Therefore, we arrive at the following non-trivial relation between hypergeometric functions: 

\ (l-e)(l+2e) 



/ 1,1 + 5,1 + 26 

' M l + e,2-e 



3^2 



r2(l+2e)r(l-£)r(l+3e) ^ f 2e,3e 



r2(i + £)r(i + 4£) 



1 



(4.36) 

The same relation can be found by considering the integral (A. 23) from Ref. ^3|, which 
they denote as 6*0(1; 1, 1; 1, 2). This integral is nothing but £4(1, 1, 1, 2, 1) in notations of 
Ref. ||3l]], which does not belong to the set of master integrals. By using recurrence relations, 
it can be reduced to the two-loop vacuum integral with equal masses ( |4.1| ) and some trivial 
part, 

1 1 

(4.37) 



£4(1,1,1,2,1) 



-lvLlll(l,l,l)- 3^3^^_;^^^_^^^ 



Ref. 54 



The result for VL111(1, 1,1) in terms of hypergeometric function is given by Eq. (33) from 

. (4.38) 



VL111(1, 1, 1) 



;i-.)(l-2.).2 



2 THl + e) 2^ ' 







Cr 





Considering the difference of the two results, one given in [^] and another presented in 
Eqs. ( [4.3?1 )- (|4.38| ), we reproduce the relation (|4.36|) . 
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It is interesting that the same function ( [4.36| ) can be obtained from Eq. ( [4.26| ), at special 
values of the indices a = 1, i^i = z/2 = 2: 



J011(l,2,2) 



1 



l,l + e,l + 2e 



'l-6){l + 2e) ' ' ' ' 



f + e,2-e 



1 
3^2 



r2(l + 2£)r(l-£)r(l + 3£) / 2e,3e 



(4.39) 
(4.40) 



where 

p'^=m^; m=l 

and we have used relation ( [4.36| ). For the same integral J011(l, 2, 2), applying Fade approx- 
imants calculated from the small momentum expansion with the PSLQ-based analysis 



(for details, see [27, ^), we restore several terms of the 5-expansion: 

J011(l,2,2) = |C2-4C3 + ^'3C4-e3{2C5 + |C2C3}+^'{fC6 + iC|} 

-e' {6C7 + 4C2C5 + 6C3C4} + {WCs + 4C5C3 + IC2CI} 
-e' {'-fC, + I2C2C7 + f C3C6 + I8C4C5 + Id} 
+e' {l^Cio + I2C3C7 + QC! + 6CIC4 + 8C2C3C5} + 0{e'). 
Using ( [4.41| ), the result for the 2-^1 function occurring in ( [4.39| ) can be deduced as: 



2e, 3e 



1\ ^ F(l + £)F(l + 4£) 
4J ~ F(l + 2£)F(l + 3£)' 



(4.41) 



(4.42) 



Therefore, we get 



1 



"3-^2 



;i-e)(l + 2e)- 

J011(l,2,2) 



+ l + 2e 




"F(l + 2£)F(l-£) 


| + e,2-e 


4; 3^2 


[ r(l + e) 



1 

3^2 



F(l + 2£)F(l-£) 

ra + e) 



- 1 



(4.43) 



(4.44) 



We note that the results similar to Eq. (|4.42|) are usually obtained using Bailey's cubic 
transformations for hypergeometric functions. 



Let us remind that for the on-shell integral JO 11 there are two master integrals [^g of 
this type, J011(l, 1, 1) and J011(l, 1,2). Using the fact that the integral ( [4.44 ) is a linear 
combination of these integrals, we obtain 



J011(l,l,2) 



F(l + 2£)F(l-£) 
'6^2(1 -2£)F(l + e) 



3e) 



J011(l,l,l) 



(4.45) 



Alternatively, one can consider other independent combinations of these master integrals: for 
example, J011(l,2,2) and [J011(l, 2, 2) + 2J011(2, 1, 2)] (see in H). While J011(l,2,2) 
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is given in Eq. ( [4.44| ), the second combination is less trivial and cannot be represented in 

ng re 

£3^2 



terms of the F- function. Using relation (see, e.g., in [^ ) 

l,l + e,l + 2£ 



3-^2 
we obtain 



2, 1 + e, 1 + 2e 
I + e, 2 - £ 



+ e,2-e 



(1 



l,l + e,l + 2e 



3-^2 



J011(l,2,2) + 2J011(2,1,2) 



1 



e{l + 2e) 



3-^2 



1, 1 + £, 1 + 2£ 



(4.46) 



(4.47) 



This is the same 3F2 function as in the case of D3, whose e-expansion is given in ( [4.32|) . The 
lowest terms of the e-expansion of Eq. ( |4.47D coincide with Eq. (9) of Ref. which has 
been obtained via PSLQ analysis, based on the small momentum expansion. 



5 Conclusions 



In this paper we have studied some important issues related to the e-expansion of Feynman 
diagrams in the framework of dimensional regularization [|l|. 

For the one- loop two-point function with different masses mi and m2, the analytic con- 
tinuation (|2.17| ) of known results |2c 



2J] to other regions of interest has been constructed. 
In particular, explicit formulae, ( |2.18D and ( ^.21[ )-( p.23[ ), relating the log-sine integrals and 
the generalized (Nielsen) polylogarithms have been obtained to an arbitrary order. 

Then, we have examined some physically-important examples of the e-expansion of the 
one-loop three-point function. For the cases of the off-shell massless triangles, as well as 
for a specific on-shell triangle diagram with two different masses, the e-expansion and the 
procedure of analytic continuation is, basically, similar to those for the two-point function. 

However, the situation for another interesting example, a massive triangle loop with 
Pi = P2 = 0, appears to be much more complicated. We have shown that all even terms of 



the expansion, i.e. the coefficients of e^', can be presented in terms of the log-sine integrals. In 
the meantime, for the odd terms (starting from e^) the result does not seem to be expressible 
in terms of known functions (related to the polylogarithms), although its one- fold integral 
representation looks reasonably simple. Nevertheless, for some special values of p| we have 
identified the corresponding contributions in terms of known transcendental numbers, 
expanding the hypergeometric functions and using the PSLQ procedure pp|. 

In particular, we have found that one extra term should be added to the odd weight-5 
basis considered in pSf. This new constant (|3.22|) can be related to a special case of the 
multiple binomial sums ( p,.ll[ ). We also constructed the even basis up to weight 5 (see 
Introduction and Appendix B.l). 

Then, we have obtained a number of new results for the higher terms of the ^-expansion 
for certain two-loop (mainly, on-shell) two-point functions (see Eqs. ([4.5|) , (|4.15|) ) and three- 
loop vacuum diagrams, Eqs. ( |4.1CI| ), ( [4.24| ), ( [4.33] ). The considered examples serve as the 

In particular, the obtained 



master integrals in the analytic computer packages [pTl 47, 48 



£-terms of three-loop vacuum diagrams are important in the four-loop-order calculations. 
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Within these calculations, a new relation ( ^.36| ) between 3F2 and 2F1 hypergeometric 
functions of argument | was found. Moreover, the analytical solution Eq. (|4.44|) makes it 
possible to express one of the on-shell integrals, J011(l,2,2), in terms of F-functions. In 
this way, we have reduced the number of non-trivial master integrals (|4.45| ) used in the 



package E? 



One of the interesting issues related to the odd and even basis construction is the con- 
nection between the generalized log-sine integrals and the multiple C, values which have 
appeared earlier in Refs. ^ , in the study of the connection between knot and quantum 
field theory. Namely, our Eq. (|A.10|) shows that 1^5,3, Ct.s and (^3,7,3 are connected with the 



combinations 



Lsr(f)-fLs,-i(f) 



the corresponding multiple C, value. 



(plus ordinary (^-functions), where j is the weight of 
Besides this, we have found that the constant U^^i (see 



in 0) can be expressed in terms of the even basis (|A.16|) . 
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A Polylogarithms and related functions 

A.l Polylogarithms and log-sine integrals 

The polylogarithm ^ij{z) is defined as 

When the argument z belongs to the unit circle in the complex z-plane, z = e^^, we get 
(see, e.g., in O]) 



Li2; (e'^) = Gl2; {9) + iCl2, {9) , U^i+i (e'^) = 012^+1 {9) + iGl2i+i {9) , (A.2) 

where Glj {9) is proportional to Bernoulli polynomial of the order j, Bj{9 j {2-n)) (see Eq. (22) 
in p. 300 of fllll), whereas Clj {9) is the Clausen function, 

n - sing / dein^'-^e ^ sin(fcg) 

^^^'^^^ - "(27^/ l-2ecos^^ + e"ti^^' ^ ^ ^ 

n - 1 / de(e-cosg) ln^'e _^ cos(fcg) 
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Below we list some useful properties of the Clausen function (for details see in p!l|): 



CI 



21+1 



\ (l - 2-^') (l - 3-^') CI2.+1 (f ) = -1(1 - 3-^') 



and 



CI 



2/ 



(1 



)l-2« 



CI 



21 



2tv 



When we consider the imaginary and real parts of Li,- ( 1 — e^^ ) , also the log-sine function 



Ls,(e) 



and the generalized log-sine function, 



d0 In^ 



/-I 



2 sin — 
2 



(A.5) 



-fc-i 



2sin- 



(A.6) 



. — e 



get involved. 
Re Li3(l-e' 
Im 113(1 
Re Li4 (1 
Im Li4 (1 
Re Lis (1 



i [C3 - CI3 (6)] + \9He, 



1 ai 



24^ 



iLsi^^(e) - i^^Ls3(^) + \eni + HCa - CI3 m k - ^^9 
-iLs4(^) + |LS3(^) k 



\eil 



Im Lisfl— e 



-^h4\e) + 1 [C5 - CI5 m + \hsf{e) k + i [Cs - CI3 {9)] I 

-i^Ls3(^) I0 + ^6'Ls4(^) - 3^6*^(3 + ^^^^0 ^ ik^"^^^' 
^Ls5(^) - 3^Lsf (^^) - iCl4 (^^) - |Ls4(0) k + iLs3(^^) ll - 



Re Lifi ( 1 — e 



^Ls5(^) - ^,Lsf\e) - \Ch (0) le - |Ls4(0) k + iLs3(^^) l] - |( 

1/3/4 I 1/3/^ /_ 1/j2t^ M^ I 1 /33?2 1 ^5 

iLs^^V)- ^Lsf V)+ 3^ [C3-CI3 (^^)] /,^+ I [C5-CI5 m k+ iLs«(e) ll 

-3^Ls('V) - |^Ls3(^) ll + 3^^Ls4(e) /e - ^^Lssl^) + ^^Lsf V) 
- ]^^^C3^e - ^^^Ls^^^(0) + ^6'^Ls3(6') - 3!^^^/^ + ^ 



Im Lig^l- 



23040 ' 



,^Ls^')(^) - y|jLs6(^) - 3^Cl6 {9) + ^Ls5(^) - ^l^^fie) I 



- j^Ls4(6') - ICI4 {0) ll + 3^Ls3(6') ll - ^Cl2 (6*) - -^Ol] + |6'C3/, 



+|^Lsi^^(e) - l,ehst\o) + i^^Cs + 4k^'Ls4(^^) - l.e^h^m I 

I 1 /1373 1 /d3a 1 /i5; 

+ 144^ - 96^ C3 - 1920^ 



(A.7) 



where /« = In 



2 sin 



— TT < 9 <t:. Results for Lij [l—e j with j < 4, as well as an outline 
how to get results for higher j's, can be found in IjlllPI 

^^A factor ^ is missing in front of hs^iO) in eq. (49) on p. 298 of |11 , whereas his Eq. (6.56) is correct. 
In eq. (7.67) of as well as in eq. (36) on p. 301, the coefficient of log^ (2sini6') CI2 {0) should be —\ 
(rather than +|). 
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According to the definition ( |A.5| ), the following integral and differential relations hold: 



LsJ-VtV) = O'^^"^) - W d0 ^'^-^Lsf (0) , ^^;:^{&) = O'^Lsfie) (A.8) 



d^ 



d^ 



In addition, using the fact that Ls2(^^) = CI2 (0) and taking into account the integration rules 
for the Clausen function, 

e e 

Chn {0) = Jd(j) Chn-l (0) , CWl (6) = C2n+1 " / d0 Chn (0) , 



one can obtain Eqs. (7.52) and (7.53) of |TT|. In particular, the function Ls^'^2(^) is always 
expressible in terms of Clausen functions, 

Ls^^)(^^) = ^^Cl2 (^^) + CI3 (^) - Cs, 

Lsf\e) = 02ci2 (0) + 20CI3 (^) - 2CI4 (^) , 

Lsf (6) = e^ch (e) + 3e''ch (o) - eecu (o) - 6CI5 (e) + 6C5, 

Ls^"^^ {6) = 24C16 {e) - 24^Cl5 (6) - ud^cu {0) + 4e^ch {0) + e^ch [e) , 



etc. 



It is also known that the values of Lsj(7r) can be expressed in terms of ^-function, for 

-i7rC2, Ls4(7r) = IttCs, Ls5(7r) : 



any j (see in [PT|): 
Ls2(7r) = 0, Ls3(7r 
etc. Moreover, 



■f7rC4, Ls6(7r) = f7rC5 + fvrC2C3 



Ls3(f) 



--C 



72' 



For 9 = ^ there are some relations among the Clausen function and log-sine integrals. 



cud) = 



(A.9) 



The result for CI4 was given in [^. Moreover, the quantities Ls^-^'^''^) can be expressed 



in terms of 7rLs,_i( ^ ) and Euler-Zagier sums. 



Ls 



5 V3 



l^Ls4(f)-fC5-iC2C 



2"bZ^35 

2 



Ls^-(fj = |.Ls5(f)+l^C6 + 2C3, 
Ls?^(f) = i^Lsed) - - TC5C2 - ^(4(3 

T „(1) f 7t\ _ l„T„/^7r\| 140879 /- I 285 /• A , 15/-2A i 21 /• 
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^vrLssf f 



192 



) /■ 945 /• A 

-<^9 — r^7<,2 



71015 a /• 12915/- A Qcz-S 



I-ttT « / 1 -U 14059623 A , 51765 A y- , 57105 /- 2 , 4305 /-2a , f^onA /• i 3075 / 



512 



,11 



26775 



C9C2 



17955 



76115403 ^ 

CsCi - 315CIC2 - if^C5,3C3 - ^C: 



16 "55 ' 4 
8875377 , 



128 



r t 

"S8<,3 



348705 



C7C4 



1917405 
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3,5,3; 



C6C5 

(A.IO) 



where Ca,fe and Ca,fe,c are defined in Eq. ( p..8|) . The result for Ls5^''(^|j was given in p8|. We 
also list some other results relevant for the weights 5 and 6: 



Ls5(f) + tvrC4, 



L4'^(f) = -^Ls6(f) + iC2Ls4(f) + |7rC5 + ^vrC2C3, 



LsfYf 



kLs5(f) + lFC6 + 4C| 



(A.ll) 



The results for Ls^^^(f) and Lsf^(f ) were given in [||]. All these results (^-(^[TT]) have 
been obtained by the rSLQ procedi 



lure. 



Using one- fold series representation ( |1.10| ) for Ls^-^"* { ^ ] and the corresponding expression 



for Lsj (f ) (see Appendix A in HQ) 
Ls.(f) = 



(2fc)! 



1 



ivr? - 1)! y 



16. 



(A.12) 



we obtain rapidly- convergent series, which can be used for the high-precision calculations of 

C5,3, C7,3 and (3,5,3- For example. 



2^-3-5 



5,3 



^ ~ {k\f 1 
26 



TT 



7 [26^tt(2A;)!A;8 
17 - 8287 - 77^ 5-19 



7I: 



1 



{2k)\ 



k=0 



{2k + iy {k\y 



16 



-C5C3 - 



(A.13) 



28.34.52.72 7 — 2-3-7" 

For the calculation of (-function with an arbitrary accuracy, an algorithm elaborated in 
Ref. 1^ can be used. 

Considering Ls^-'''(^|^ and Lsj''(7r) (j = 4, 5; / = 1,2) we find that most of them (except 
which is one of the ewen-basis elements) are connected with Lij , 



for Lss^"* 
Ls«(f) 

Lsl'V) 



-I In^ 2 + In^ 2 - iC3 In 2 + f C4 + ivrLsg (f ) - |Li4 (i) , 



i,ln^2 + AC2ln3 2 



16 



i2|C3ln^2 



15 



Li4(i) ln2-|C2C3 



Ls 



(1), 



77 



+ ivrLs4(f ) - WC5 
-| In^ 2 + C2 In' 2 - IC3 In 2 + f (4 - 4Li4 (f) , 

-I In^ 2 + 2C2 In^ 2 - f (3 In' 2 - 12Li4 (i) In 2 + §(2(3 + §(5 - 12Li5 (i 



Lsf (tt) 



-ivr In^ 2 + 277(2 In' 2 - 777(3 In 2 + f 77(4 - 87rLi4 (} 



(A.14) 



general factor, (—1)", in definition of Ls„(0) is missing in Appendix A in 
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All relations ( |A.9| )-( |A.14| ) have been obtained using the PSLQ procedure [^0. We give 
some relevant numerical values, 



Ls3(f 



- -2.03357650607205460091206896970 . . . , 
^ 6.003109556529006567309305614033 . . . , 

- 0.9889445517411053361084226332284... , 

- -24.0143377201598359235946799181446. . 
^ -0.1268132428355886971002322996611 . . . 

As an example of application of the even basis, we consider 

Ua,b = C{a,b] -1, -1) . 



Ls5(f 



(A.15) 



The lowest basis element of alternating Euler sums [§] which cannot be expressed in terms 
of (j, In 2 and Lij (^|^ or their product^ is f/5,1. This constant appears in the e-expansion of 
three-loop integral B4 (see details in ||68|) which is connected with the charge renormalization 
in QED (Eq. (51) in as well as in the 0{1/N^) contribution to the critical exponent in 
the large- limit in (3 — 2e) dimensions (see Eq. (23) in |69]). The PSLQ-analysis, with an 
accuracy of 300 decimals, yields the following relation: 



39 ^^6 {2 



^vrLssff) + i^Liefi) + i^Lisfi) In 2 + ^Li^fi) 111^2 



13 



13 



13 



+ lis 2 - ^Ce - M + iC5 In 2 + ^(3 In^ 2 - ^(2 In^ 2 . (A. 16) 
For the procedure of the analytic continuation (in Section 2.2) we need the generalized 



(Nielsen) polylogarithms which are defined as (see, e.g., in Ref. [35]) 



Sa,b{z) 



(-1) 



a+b-1 



ln"-^ein'(l-^0 



(a - 1)! 6! 

In particular, 

Sa,i{z) = Li 
The following integration formula is useful: 

Sa,b{z) - i-lfSaA^/z) 



I 



^a+l{z) , Z — Sa^biz) — Sa~l^b{z), 



dz 



So,b{z) 



(A.17) 



-1)' 



6! 



ln''(l - z) 



(A. 18) 



fc=0 



{p-k)\ 



Sa+l+k,b{^) + i-lf+'Sa+l+kA'^/z) 



(A. 19) 



21r 



The result for Ls4^''(7r) coincides with Eq. (7.71) of whereas in Eq. (7.145) for Ls5^''(7r) the term 
+27rC2 In^ 2 is missing, and the term +377r^/360 should read 37r^/40. There is also an error in Eq. (7.144) of 
Iprf: Ls^^^(27r) should be equal to -137r5/45, rather then to 77rV30. 
22Two other constants are C(5, 1, -1) and C(3, 3,-1), sec Eq. (17) in |]. 
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Starting from Eq. ( |2.18| ) and applying relations ( |A.8| ) and ( |A.19| ), we arrive at 



(1) 



Ls j(9) — Lsj(7r) 



2^J(J + 1) 



\riP(-z) ^'^'^ 



p=0 ^ f k=l 



+2(-iy(j - 1)! ^ ^A:5,+,,_,(-l) 



fc=i 



2fc 



(A.20) 



Substituting 6* — vr — icrln(— 2;) (see Eq. ( p.l5|) ) and using Eq. ( p.l8| ) we can rewrite this 
result in a different form, 



LsYiW-LsYi(vr) = vr 



Ls,(^^)-Ls,(7r) 



+ 



ln^+^(-^) l-(-l) 



2^(j + l) 



1)* 



+2(-iy(j - 1)! 5: ^kS,+,^,.,i~l) 



k=l 



p=i ^ fc=i 



A. 2 Auxiliary function 

It is useful to define another generalization of the log-sine integrals as 



(A.21) 



Lscij(^) = -J d(j)W-^ 



2 sin ■ 



In^-^ 



2 cos ■ 



(A.22) 



When z = 1 or j = 1, it reduces to ordinary Ls functions, 

LsCj_i(^^) = Lsj(6') , Lscij(^^) = — Lsj(7r — 9) + Lsj(7r) 
There exists an obvious symmetry property, 

LsCjj(6') = — LsCj_j(7r — 9) + LsCj^j(7r) , 



(A.23) 



(A.24) 



which implies that the functions LsCj ,,(^^) with i > j are not independent. The values of 
LsCjj(7r) can be extracted from Lewin's book |jTl|, Eqs. (7.114)-(7.118). We present here 
some of them 



Lsc2,2(7r) = i7rC2, Lsc2,3(vr) = -ivrCa, Lsc3,3(7r) = -^71(4, Lsc2,4(vr) = §7rC4 • (A.25) 
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Using sin = 2 sin | cos |, the following duplication formula can be derived: 



2 2 ' 

k-l 



Ls,(2^) = 2{k - 1)! J2 ^TTtA ^Lsc,+i,fe_,(0) . (A.26) 



i\{k — 1 — i 

For special cases k = 2,3,4,5 it gives, respectively, 



Ls2(^) -Ls2(7r-^^) = iLs2(2^), 

2Lsc2,2(^) = iLs3(2^) - Ls3(^^) + Ls3(7r -9)- Lss{n) , 

3 [Lsc2,3(^) - Lsc2,3(7r - e)] = |Ls4(2e) - Ls4(^) + Ls4(7r -9)- ^Lsiin) , 

6Lsc3,3(^) + 4 [Lsc2,4(^) - Lsc2,4(7r - 9)] = iLs5(2^) - Ls5(^) + Ls5(7r -9) + f 71(4 • (A.27) 

We can see that for odd values of k one would always obtain representations of the functions 
Lsc(fc_|_i)/2, {k+i)/2{0) in terms of the Lsc functions with i 7^ j. This means that the functions 
Lscjj(^) with i = j are not independent. Therefore, it is enough to consider only the 
functions with i < j. In particular, up to the level k = 5 only two new functions are needed, 
in addition to ordinary Lsj(6'): Lsc2,3(6') and Lsc2,4(6'). 
For a particular point 9 = ^, Eqs. ( |A.27| ) yield 

Lsc2,2(f) = -|Ls3(7r) = |7rC2 , 

Lsc2,3(7r) = -iLs4(7r) = -|7rC3 , 

Lsc3,3(f) = -|Lsc2,4(7r)- i^Ls5(7r) = -^7rC4, (A.28) 
whereas for 6* = | we get 

Lsc2,2(f) = |Ls3(f ) - iLs3(f) + , 

3 [Lsc2,3(f ) - Lsc2,3(f )] = |Ls4(f ) - Ls4(f) - |7rC3 , 

6Lsc3,3(f ) + 4 [Lsc2,4(f) - Lsc2,4(f )] = iLs5(f ) - Lssd ) + f ""^^ • ^^■'^^^ 
For these particular values of 9 the PSLQ procedure yields 

LsC2,3(f) = -K3-^Ls4(f) + iLs4(f), 
LsC2,3(f) = -^7rC3 + ^Ls4(f), 

L-^,4(f) = l<4 + !C2Ls3(f)-i^Lsl^)(f)-^Ls5(f) + |Ls5(f)+iLsrO 
^Mt) = li-C4 + ¥C2Ls3(f)-i-Ls«(f)-iLs5(f) + iLsf(f 
Lsc3,3(f) = -l-C4 + iGLs3(f)-|vrLs«(f)-iLs,(f)+iLsr(f 
Lsc2,3(f) = |7rC3-|Ls4(f)+ 2014(1). (A.30) 
These results obey the conditions ( [A.28|) -( |A.29| ). However, no relation has been found for 
Lsc2,4 (f ) ^ 0.30945326106363459854315773429895417071 . . . , (A.31) 
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which represents a new transcendental constant appearing at the weight-5 level of the even 
basis. 

If we introduce a variable z = e'"^, we see that 

In (2 sin I) ^ ln(l -2) - iln(-^), In (2 cos f) ^ ln(l + ^) - | Inz, 

Therefore, the analytic continuation of Lscj^j (6) is related to the integrals of the type 

dz 



y^ln°(l + 2)ln^(l-z) In^ 



z . 



This is nothing but a particular case of the harmonic polylogarithms |7^. It is known 
that some of such functions of weight 4 cannot be expressed in terms of polylogarithms 
(or Nielsen polylogarithms) of specific types of arguments, like ±z and |(1 ± z) (see in 
Ref. |]T2[). For example, in our case the analytic continuation of Lsc2,3(6') involves an integral 
J ^\n^{l + z)\n{l- z). 



B Hypergeometric functions 
B.l Procedure of the e-expansion 

Let us briefly describe some technical details of obtaining terms of the e-expansion of hyper- 
geometric functions 2-^1, 3-^2 and 4F3 given below. All of them belong to the type 



Ai + aiE, . . . , Ap + apE 

B + \ + hE,Ci + CiE, Cq_i + Cq_iE 



(B.l) 



where Aj = B = 1, whereas Cj take the values 1 or 2. More generally, we shall assume that 
Aj, B and Cj are positive and integer. This function ( [B.!] ) can be presented as 



1 {Ai + aiE)j . . .{Ap + apE)j (2i3 + 1 + 2fe) 



j! (Ci + cie)j . . . (Cq_i + CQ_i£), (B+l + bE), 



2j 



(B.2) 



where {a)j = r{a + j)/r{a) is the Pochhammer symbol, and we have used the duplication 
formula (2/3)2, = 4^(/5)j(/5 + |),. 

To perform the e-expansion we use the well-known representation 



[1 + aE)j = j\ exp 



°° ^-aE)'' 



E 

k=l 



SkU) 



(B.3) 



which, for ^ > 1, yields 



{A + aE)j = {A)j exp <^ - 

I k=l 



-aE) 



[Sk{A + j 



Sk{A-l)] 



(B.4) 



where Sk{j) = Z]z=i ^ is the harmonic sum satisfying the relation Sk{j) = Sk{j — I) + j 



-k 
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In Ref. it was demonstrated that the multiple binomial sums ( |1 . 1 1| ) , or their certain 
linear combinations, can be expressed in terms of transcendental constants elaborated in 
28[| . The only values of A, B and C for which each term of the e-expansion of ( |B.2| ) can be 



expressed in terms of the multiple binomial sums (|1 . 1 1|) are B = 1, A = 1 or 2, C = 1 or 2. 
For other values of the indices binomial sums of a different type appear, which are connected 
with series like 

|i;p;r^n;^.("-i)ris.(2«-i)f, 

where p and q are positive integers. In such series, new transcendental constants appear. 
For example, for p = 1, q = 2, i = j = we get (see in W^) 



^J2n)! (2n + 1)2 3 3 V v 

where G is the Catalan's constant. 

The values of multiple binomial sums ( p.2|) or their combinations up to weight 4 are 
presented in Appendix B of Ref. Most of the sums can be expressed in terms of 

the basis elements. There are, however, a few weight-3 and weight-4 sums which are not 
separately expressible, but only certain linear combinations of them, namely: 

+ ^-;2;l5 ^-;1;2 + ^-;2;2) ^1;1;1 + ^li2;l ) "^-ll;! ~^ ^^~\l,2;l ~^ ■ (B-5) 

The situation at the weight-5 level is similar. Again, most of the sums can separately be 
expressed in terms of the basis elements, except for a few ones, which are expressible only 
in certain linear combinations. There are 4 linear combinations independent of xs ( |3.22| ), 

^l;l;l "T ^l;2;l) ^2;1;1 + ^2;2;1 5 ^1;1;1 "I" <J^1;1,2;1 + ^^l;3;l5 

+ 3e:;|i + gejXi + 8s::l;^;i + 6s:g^i, (b.6) 

and 3 combinations which involve Xs? 

^-;l;3 + ^-;2;35 ^-;1;2 + '^^-\3;2 + '^^-\i,2;2j ^1;1;2 + ^l|2;2- (^-7) 



Using relation (|1.9| ) and the results of Ref. [^], is it natural to expect that the sums 



not only with k = 1 but also with k = 3, are connected with the odd basis, whereas the 
sums (|1.11|) with k = 2 and /c = 4 should be associated with the even basis. 

To check this conjecture, we have performed PSLQ-analysis of the sums ( |1.11| ) with k = 2 



up to weight 50. At the weight-5 level, one needs to add two new elements to the even basis 
of One of them can be associated with Lsc2,4(-|) (see Eq. ( A.31| ) in Appendix A. 2) 



whereas the second one, xs; is given in Eq. ( |3.28| ). In this way we restore, accidentally, 
the "empirical" relation Nj = 2^ for the j = 5 level of the even basis. Then, we find that 
the expressions for the sums ( |1 . 1 ID with k = 2, contain all elements of the even basi^. 



23n 



The values of some sums ( 1 . 1 1| ) with fc = 4 up to weight 3 can be extracted from the resuhs of (TI] 



This is connected with the fact that both Ls^^^(^) and Lsj^''(7r), for j — 4,5, are expressible in terms 



3 

of the same constants (see Appendix A.l 
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However, the question about possibility to express all elements of the even basis in terms 
of only binomial sums (|1 . 1 1|) is still open. As in the k = 1 case, most of the k = 2 sums 
are separately expressible in terms of the even basis elements, except for the same set of 
sums, which can be expressed only in certain linear combinations^: namely, for the weights 
3 and 4 the combinations given in Eq. ( [B.5| ) and at the weight 5 the familiar combinations 
( p.6|) and ( |B.7| ), everything with k = 2. Worth noting is that the combinations ( p.6|) do 
not contain X5j whereas the sums ( |B.7| ) do not involve Lsc2,4(^|). In total, we have found 
8 elements which involve X5'- besides the combinations given in ( p.7|) (with k = 2) and the 
sums S^:}.2(2), S}:}.3(2), Til'.~.-^{2) (which have already occurred in the k = 1 case), there are 
two other sums, Zl}c.4(2) and SI;J.4(2) (for k = 1 these two sums are expressible in terms 
of Cs, C2C3 and 7rLs4(f)). 

We would like to mention some issues related to the high-precision calculation]^ of the 
elements of our basis, needed for the PSLQ-analysis. Most of the constants of the basis 
up to weight 5 can be produced, with desirable accuracy of 1000 or more digits, using the 



algorithm described in Appendix A of Ref. [33|. The Xs aiid Xs constants are originally 



defined in terms of rapidly-convergent series, see Eqs. ( |3.22| ) and ( |3.28| ). Therefore, only 
three constants, L sf^(^), Lsc2,4(f) and Ls^^^(f), may require additional consideration. In 
principle, the corresponding binomial sums may be used for precise numerical calculation. 
However, specifically for Lsf^ the following representation in terms of inverse tangent 



integral (see in ||rT[) happens to be very useful: 



Ti5(;^) - E 3, i ■ + 1)5 = 2k^ 3 + In^ 3 + ^vrCs In 3 + igvrC4 



J- 
5 

'288 



1^0 3^' (2j 

Ls2(f ) ln=^ 3 + ^Ls3(f ) In^ 3 - ^,Ls,{^) In 3 - l,Ls,{^) In 3 
-^C2Ls3(f ) + ^Ls5(f) + iLs5(f ) + ivrLs«(f ) - ^Lsf (f ) .(B.8) 

For numerical calculation of Ls^''\6) the following algorithm is useful. Substituting y = 
sin|, we get from the original integral representation ( |A.(j| ) 



sin(0/2) 



Lsf\0) = -2'+' I ^L=(arcsiny)'=ln"-'=-i(2i/). 



Expanding 



{2r)\y 



Such interrelation is not trivial. For example, for fc = 3 the number of such combinations of sums is 
reduced, so that the first irreducible combination arises only at the weight 4. Besides this, the element X5 



is not reproduced by series like (3.22) with A; = 3. 

^^AU such computations have been performed with the help of David Bailey's MPFUN routines. 
The latest version of these codes, including the documentation and original papers, is available from 



tittp : //www . nersc . gov/ ~dhb/mpdist/mpdist . html 
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and (arcsinj/)'^ (see below) in y and using 



x^Iti^x dx = {-if b\ x^-^^Y. 



Inx 



.b~p 



-o(6-p)!(a + l)P+i ' 

we obtain the multiple series representation suitable for numerical eveluation. 

The expansion of (arcsiny)^ is discussed in Ref. |7^. The generating expression is 
exp(a arcsin?/) = J2'^o{bpy^ / p\) , where bo = 1, bi = a and 



62, = a\a' + 2'){a' + A^)---{a^ + i2k-2Y), k>l. 

Expanding exp (a arcsin y) in a power series in a and equating coefficients of on both sides, 
the Taylor series for (arcsin y)'^ can be deduced. For example, 

1 



2! 
1 
3! 



(arcsin y)" 



[a-Tcsiny)' 



.^0 (2r + 2)! ' 
-/ I 1 \ (2r)! 

2^1 +32+ ^ (2r- 1)2^ (r!)2 4'-(2r + 1) ' 



— (arcsin y)^ 



r=l \ 

^ / 1 1 



1 



4'"(r!)2y2r+2 

(2r)2 y (2r + 2)! 



B.2 The 2F1 function 

Here we present some relevant results for the hypergeometric function of the type 



/ 1 + aiE, 1 + 
'^'\ l + be 



(B.9) 



Especially we are interested in its expansion in e, which basically corresponds to the calcu- 
lation of the derivatives of 2^1 function with respect to the parameters. 

There are some special cases when the result is known exactly. For instance, when 
Oi = —02 = a and 6 = we have (see Eq. (94) in p. 469 of p5||). 



„ / 1 + ae,l — ae 



sin^^ 



sin {2ae0) 
ae sin(2^^) 



(B.IO) 



When ai = b = Q {a2 = a) we get (see Eqs. (^)-(^)) 



2F,^ I 



sin^^ 



2l+2ae sin 0fcOS 0)1-2'^= 



j=0 



(B.ll) 
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i.e., we know all orders of the e-expansion, in terms of the log-sine integrals ( |A.5| ). For 



= |) we get Eq. ([4.16|) . Another interesting case is ai = 0, 02 = 2& = a, for which we find 



/ 1,1 + ae 



1 + ae 



(B.12) 



Moreover, in a more general case, when ai = while 02 = a and b are arbitrary, the 2-F1 
function ( [B.9| ) can be reduced to an incomplete Beta-function (see in [0). As a result, we 
get the integral representation 



/ 1,1 + ae 



l + 2be 



smt^ cos 6* 



d0(sin0)2^"(cos 



\2ae~2he 



(B.13) 



The ^-expansion can be written in terms of Lsc function (see Appendix A. 2). 

There are some further results available for the case when z = \ {Q - 
a\ = h= a and 02 = fa we get 



/ 1 + ae, 1 + lae 
l + ae 



2ti 



r(2 + 2ae) 



gf (i+ae) r(l + ae)V^{l + \ae) 



f). When 



(B.14) 



Using the procedure described in Appendix B.l, we obtain a few terms of the e-expansion 
for general values of Oi, 02 and h: 



/ 1 + aiE, 1 + 

2-ri 3 



+ he 



l\ 2(1 + 265) 



X ■ 



ivr + ei2A, - 56) J-ff^ib - A,y [Ls,+2(f ) - Ls,+2(vr) 



+e 



|(A2-A2)(2Ai-56)C2Ls2(f 



lb' + ±{Al-A,)-f-^bA,\ 7^(2 + e' 
+ {W2MAI - A,) - 1|6^ - f|l6A? + 116^2 + fb^A,) <3 
+ (-f (2Ai - 5b){Aj - A2) + ^6(6^2 + 96Ai - 426^)) Ls4(f 



+e'[^{Al - A2){2A, - 56)(llAi - 356)C3Ls2(f ) - ^(A^ - A2){2A^ - 5b)\ [Ls2(f 

+1(6 - Ai) (-f (2Ai - 56)(A2 - A2) + 66A2 + 96^1 - 426^) C2Ls3(f ) 

-lib - Ai) (-8(2Ai - 56) (A? - A2) + 66A2 + 962^1 - 426^) [nLs[^\^) - flsf (f 

_l_f 2091,4 17531,3 /1 , 10399 ^2 4 11069 ^2 42 , 10061 ?, 4 3 10799/, /i /i 
+ 1^^'' ~ W'' ^1 + "648"'' ^2 - 7296"'' ^1 + "648"''^1 ~ "648"''^1^2 

,305424 5491 44 I 1 42\ /■ I 11 1,4 53x,3 4 , 2937,24 433l242 
+ ^A^A2 - 1296^1 + 1296^2 ^^Qi + [y'^ " is'' ^1 + ^'' ^2 " Yga'' ^1 



+f>6A? - II6AA2 - f (A? - A2)A? Ls5(f ) 



+ 0(5' 



(B.15) 
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where Aj 



=1 1 



and we take into account that there are some relations among Aj, 



2^3 + Al- 







2Ai + A\- 2AIA2 - Al = Q . 



At the weight-5 level of the odd basis, we have observed that the constants '7rLs4^'' ^— 



and Lsf (f 



vrLsi'Y^) -!Lsr(^ 



^4 IT 



5 V 3 



Vl3 

Moreover, this 



appear only in the combination 

combination arises only in the sums like Sa4 i- It is natural to expect that similar "junctions'' 
also happen for higher generalized log-sine integrals. 

We also mention some useful results for the contiguous 2-^1 functions. One of them 
Eq. ( [4.42| ), is rather essential for the results of this paper. Another interesting case. 



2F1 



|-2e 



2'-^^ r(|)r(| 



2e) 



31-3- r(i)r(| 



2e) 



(B.16) 



can be obtained by a simple transformation of Eq. (114) on p. 461 of E5| 



B.3 The 3F2 function 

Here we present some results for the hypergeometric function of the type 



11 + aiE, 1 + a2e, 1 + a^e 
l + be,C + ce 



(B.17) 



with C = 1 or C = 2. 

There are a few special cases, when the considered function ( |B.17| ) can be expressed in 
terms of 2F1 functions. First of all, when C = 1 and c is equal to one of the (for example 
c = as), we obtain the function ( p.9|) . Another reduction case is C = 2 and c = 03 = (see 
Eq. (7) in p. 497 of and Eq. ( ^ ) of this paper). 

There is also a less trivial relation, valid for the case when C = 2 and c = 2b = 2a^ = 
ai + a2 (see, e.g., Eq. (22) on p. 498 of Q), 



^ H i + |(ai + 02)^,2 + (ai + 02)^ 



2F1 



1 + \aie, 1 + \a2e 
I + \iai + a2)e 



B.18) 



which reduces the 3F2 function to a square of the 2F1 function of the type ( [B.9|) . 
Another interesting relation (which follows from Eq. (20) on p. 498 of |H5|) reads 



3-r^2 



1 + flie, 1 + 02^, 1 + \{ai + a2)e 
I + \{ai + a2)e, 1 + (ai + a2)e 



2F1 



1 + ^aie, 1 + ^a2e 
I + |(ai + a2)£ 



1 



aia2S z 



2 [1 + (ai + 02)5] 



^ / 1 + \aie, 1 + \a2e, 1 
^ M I + i(ai + 02)5, 2 



[B.19) 



This relation gives a possibility to cross-check results for the e-expansion of both 3F2 functions 
presented below, together with the result for the 2F1 function. 
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We note that the results for some other values of the parameters can be obtained using 
relations between contiguous 3F2 functions, like (see Eq. (26) on p. 440 of [|^) 



^^ + ^^^)^^^[ l + be,2 + ce 



(l + ce) 3F2 



l + ai6, l + a2e, l + a^e 
^ + be,l + ce 



+ (03 - c)e 3F2 



l + ai6, l + a2e, l + a^e 
^ + be,2 + ce 



.(B.20) 



This connection makes it possible to use, as a check, another special case (see case (iii) in 
Table I of pl). 



p I 1 + aie, 1 + 026, 2 - 



gaie-f 



r (I + \{ai + a2)e) T (2 + ^(01-02)5 



X < 



23r(| + |ai£ 



3a2£:(l — a2£:)r(l + aie) 
r(l + iaie) 



> . 



r (I + i(ai-3a2)e) T (| + l{ai + 3a2)e) T (l + |(ai-3a2)e) T (| + |(ai + 3a2)e) 

(B.21) 

Again, using the procedure described in Appendix A.l, we obtain a few terms of the 
^-expansion for general values of Oj, b and c. For the case C = 2 we find 



3-^2 



1 + aiE, 1 + a2£, 1 + a^e 
^ + be,2 + c€ 



2(l + 2be){l + ce) 



x\k2 + e ICs (11^1 - 356 + c) - |7rLs2 U) {2Ai - 5b + c 



i7rLs3(f ) {2Ai - 56 + c)(Ai - 6 - c) - I [Ls2(f )]^ (2Ai - 56 + c)(2Ai - 56 - 2c) 



-L/- i 1707,4 139^4 -U 377l„ I 160l2 211 2 ,108542 5 4 

+U [ — ^bAi - ^cAi + -^bc + — 6 - YogC + -216^1 - 216^2 



^X5(2Ai- 56 + c){2Ai- 56 - 2c) (Ai- 6 - c) - |7rLs4(f ) 6(2Ai- 56 + c)(Ai - 6 - c) 



-T|^7rC2Ls2 (f ) (2Ai - 56 + c){109b'^ - 23bAi + 236c + 2c^ - 2^2 



+C2C3(fifeA? - H fecA, - f 6M, + icAj + ^,A,A2 - I6A2 + ^CA2 

_ 13 24 _ 773 L 2 , 833 l2 , 3361 l3 _ _53. J 
^27 3 81 ^1 fi4s"^"^ ^ ins"^ ^~ fiA« ■^9/1"^ 



81 



648" 



+C5 



1903 ^2 
72 



6Mi - ^6cAi - fi6A^ + icA^ + i(cA2 + ^AiA2 + iA, 



108" 



529/ 



648 



47 



324 



17. 



53 



58 



I 164 24 I27IL4 4217l„2 I 1445l2 22907 l3 , 103 „3 



162 



648 



108 



+vrLs4(f ) l6cA, - §6% + li6A? 



648 



20 ^2 



324 



243-c'^l + M|&^2 



243 



cAo 



10 
243 



A1A2 - IcAl + §^bc^ 



29 x,2. 



ilj^C + ^6^ 
81 " 486 



^C^ 
243 



+ 0(£' 



(B.22) 
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where Aj = J^^^iaj, while xs is defined in Eq. ( p.22| ). In particular, we have checked that 
this result obeys the condition ( p.l8|) , where the e-expansion of the resulting 2-^1 function 
is given in Eq. ( B.15| )). 

For another case of interest, C = 1, we find 



3-^2 



1 + aie, 1 + a2e, 1 + ass 



2 



+ beA + ce 



2(1 + 2be) 

-Ai€—be—c£ 



32 



5b 



b-cy [Ls,+2(f) -Ls,+2(vr) 



+e^7i(247 \aI - + (c - Ai)(156 + 2c) + 576' 



_475l 2 , 409 L 4 _ 13 l2 _ 17l3 , 49_ a _ 49 3 _ 16 4 
81 ''^ + 162''"^2 gOC ^^0 + ^^^CA2 gl^ 27^^ 



162 



27' 



162" 



81' 



+C2Ls2(f) (^A^bc+lAi^ + lAiA2-lAlb-\Alc+lAl-^p^ + lbA2-\cA2 + 
+Ls4(f ) (1^6' - f^Ai6c - ^Aic' + If + '^Alb + f A'c - f,Al 



9^ 9 



+^bc' - ^bA2 



+e 



\b'c - §b' - icA2 + §c' + lA, 



81 



81 



C2Ls3(f ) (Ai - 6 - c) {aaI - fAlb - 12Alc + fA^bc - SA^b^ 
-lA,A2 + f Aic' - f fee' + 36^0 + 1463 + fcAa + f - f A3 - fc^) 
+C3Ls2(f ) (|Aic3 - f A1A3 - f^Aifec' - §A,A,b - ^-§A,b'c - -^A.A.c 



+^A{bc + ^^252 + ^A{c' + A 



19 /l2^2 



70 
27 



Alb - lA\c + \A{ + f 6CA2 



80^ 3 



I „^ / 78145 4 /,„2 39377 /I 4 /, i 11069 4 /,2 1753 4 /,3 1843 4 4 „ , 11021 4 „3 



38971 



5279^^^ 
.6125 44 I 39377^^^^ _ 37901 



432 ^"i^^3 seF'^i"'^ 



11069 42^2 _ 6931 42„2 _ 1619 42 4 , 38971 /i3l , 9475 43„ 
1296 "^l" 288 288 ^l^'^^ + 2592 ^l" 648 "^1*^ 



2592 

20£ 
36 



2592 



L d _ l^/)4^ 
1296 1296 -5 



31867^2 2 , 10399 12 a . 1753 u3 „ 



1296 



I 209^4 _ 2437 4 i 3665 2 4 1223 4 , 5491 4 2 305 4 
~r Qft " ionfi''^3 ~r" 2592 2 432 ^ "i" 1296 2 .'^fi 4 



+7r 



1296 
2 



432 



36 



Ls2(f )] {§A^bc^ + ^AM + §A,b'c + ^A^A^c - f^ic^ + §A^As - fAjbc 



-§Alb' - l^Ajc' - §AlA2 + §Alb + |A?c + - §bcA2 + §bc' 



-§bA, - §b\^ + §b^A2 - §,cA, + ^c-A, - ^c^ + 4,Al - 4,A, 
+Ls5(f ) {^A,bc' - '-i^AM + f Aifo'c - §A,b' - i|AiA2C + ^^A,c' + ^A^A 

1409 A-iu^ 433 A^l.^ 2093^2^2 _ 665 ^2^^ _^ ^A% + —A^C- —A^ + — 6cA 

53 l3 



108 
162 



■Albc- i|A?6^ 



43342^2 

162 1" 324 
10191,2 2 I 293^2 



09 43^ , 109 43 583 44 , 1231;^ 4 1249?,^: 
!4"^l"+ 27 "^1*^ 972^^1+ 324 "'^"^2 162"'^ 

711,4 1019 l,:^„^ I 293 1,:^ 4 I 537,3 „ I 117,4 107„4 _|_ 95„2 4 95 „4 , 104 42 208 4 

-—bAs - -j^b c + ^b A2 + j^b c+ ^b - jq^cA^ + ^iC A2- j^c + - ^A. 



ttLsI') (f 



3t„(2)/2^ 

- 4^S5 y— 



{A^-b- c) (f|A? - f A'fe - f A'c + f Ai6c - fAife' 
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-^A^A2 + f Aic^ - fbc" + p^c + 7b^ + ^cA2 + fbA2 - f A3 - |c=^) + 0{e'')^ , (B.23) 



B.4 The 4F3 function 

Finally, we present some results for the hypergeometric function 



11 + aiE, 1 + a2e, 1 + a-^e, 1 + a^e 
^ M l + be,2 + cie,2 + C2e 



(B.24) 



There is an interesting special case, ci = ai and C2 = 02, when this function (p.24| ) 
reduces to a combination of two 3F2 functions of the type ( |B.17| ) with C = 2 (see, e.g.. 



Eq. (4) in p. 497 of [||]) 



/ ^ / 1 + aiE, 1 + a2e, 1 + ase, 1 + a^e 

e{a2-a,),Fsy | + fe, 2 + a^e, 2 + a2e 



, , / l+aie,l + a3e,l+a4e 



l + be,2 + a2e 



z IB.25) 



As in the previous cases, using the procedure described in Appendix A.l, we find a few 
terms of the e-expansion for the general values of Oj, b and Cj, 



/ 1 + a\e, 1 + 026, 1 + 03^^, 1 + (^aB 

4-^3 I 3 



f + fee, 2 + cie, 2 + C2e 



2(1 + 26£)(1+Ci£)(l + C2£) |7rLs2(f ) - |C: 



+£ 



+£^ 



LS2 



(2Ai - 56 - 2Ci) + C4 6 - If Ai + 7Ci 



-Ls3(f) (6-Ai + Ci) + 

-X5|(2Ai - 56 - 2Ci)(Ai - 6 - Ci) + 7rLs4(f ) 6(^1 " ^ " C^i) 
C2Ls2(f ) ^ (10952 - 236^1 + 236Ci - 2^2 + 2^2 



I A /• ( 1063/,^ I 1 4 4791,4 -U 1 13 42 I4 , 197/,2 

I A ( 1103 L 4 127^ 4 , 53 ^2 71^ 29 4 47 42 , 2395 8599 ^2 



+7rLs4 (1) (ifCi^i ~ ^^"^l ~ fH^^l + 1^1 + ^^2 + y||&^ + ^C2 - 
where A^ = Y^i^i aj, Cj = c-, while X5 is defined in Eq. (|3.22| ). 



0(£3) ,(B.26) 
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